PHOTOVISCOELASTICITY 


DISTRIBUTION  STATEMENT  A 

Approved  for  Public  Release 
Distribution  Unlimited 


by  Ellis  Harold  Dill  and  Charless  Fowlkes 


Prepared  under  Grant  No.  Ns( 


>f) 

j-401  by 


&  UNIVERSITY  OF  WASHINGTON 


^  Seattle,  Wash. 
for 


20020319  177 


-  WASHINGTON,  D.  C.  -  MAY  1966 


NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 


NASA  CR-444 


PHOTOVISCOELASTICITY 
By  Ellis  Harold  Dill  and  Charless  Fowlkes 


Distribution  of  this  report  is  provided  in  the  interest  of 
information  exchange.  Responsibility  for  the  contents 
resides  in  the  author  or  organization  that  prepared  it. 


Prepared  under  Grant  No.  NsG-401  by 
UNIVERSITY  OF  WASHINGTON 
Seattle,  Wash. 

for 

NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 


For  sale  by  the  Clearinghouse  for  Federal  Scientific  and  Technical  Information 
Springfield,  Virginia  22151  —  Price  $1.45 


ABSTRACT 

report  provides  a  basis  for  the  description  of  the 
mechanical  and  optical  properties  of  biref ringent ,  visco¬ 
elastic  materials,  and  for  the  application  of  these  materials 
in  experimental  stress  analysis  by  photomechanics.  The  ex- 
periment al  pro cedure  applied  to  two  plane  stress  problems. 

Preparation  of  the  model,  characterization  of  the  material, 
and  the  test  apparatus  are  also  described 
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1. 


INTRODUCTION, 


Many  problems  in  structural  engineering  related  to  the 
structural  integrity  of  solid  propellants  or  the  deformation 
of  heated  metal  and  ceramic  structures  can  be  reduced  to  the 
stress  analysis  of  a  linear  viscoelastic  material.  The 
stress  analysis  of  bodies  with  irregular  boundaries  is  such 
a  complicated  mathematical  problem  that  it  is  frequently 
more  simple  to  determine  the  stress  state  by  experiments 
upon  a  model.  This  report  deals  with  the  extension  of 
the  method  of  experimental  stress  analysis  known  as  photo¬ 
elasticity  to  viscoelastic  materials. 

The  basic  idea  is  to  manufacture  a  model  from  a  trans¬ 
parent,  birefringent  material  and  to  determine  the  stress 
state  of  the  model  by  observation  of  the  fringe  pattern  which 
occurs  when  the  model  is  placed  between  polaroids.  In  order 
to  correlate  the  stress  in  the  prototype  with  that  of  the 
model,  the  mechanical  properties  of  the  model  and  the  proto¬ 
type  must  be  similar.  Thus, the  material  of  the  model  must, 
generally,  be  viscoelastic.  Birefringent  materials  having 
the  desired  mechanical  properties  can  be  manufactured  from 
high  polymers. 

The  fringe  pattern  observed  when  birefringent  materials 
are  placed  between  polaroids  is  due  to  the  dielectric 
properties  of  the  material.  The  isoclinic  and  isochromatic 
fringe  patterns  give  information  about  the  dielectric 
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properties  at  any  instant. 

The  dielectric  properties  of  a  bi r ef ringent ,  visco¬ 
elastic  material  at  any  instant  depend  upon  the  entire 
history  of  deformation.  The  fringe  pattern  changes  with 
time,  even  when  the  stress  or  strain  is  constant. 

The  relation  between  dielectric  properties  and  mechani¬ 
cal  state  can  be  inverted.  Thus,  the  strain  or  stress  at 
any  instant  depends  on  the  entire  history  of  dielectric 
properties.  This  means  that  it  is  necessary  to  have  a 
complete  time  record  of  the  history  of  isochromatic  and 
isoclinic  fringe  patterns  in  order  to  gain  information 
about  the  present  value  of  the  stress.  This  is  the  main 
difference  between  photoelastic  and  photo vis  coelastic  pro¬ 
cedures  . 

The  applications  described  in  this  report  will  be 
limited  to  quasi-static  problems.  It  will  be  seen  that 
photoviscoelasticity  is  a  practical  tool  for  engineering 
design. 
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2. 


VISCOELASTICITY. 


2.1.  Basic  Equations. 

In  the  following  sections,  some  results  from  the 
classical  theory  of  viscoelasticity  are  presented.  The 
equations  are  limited  to  small  displacements  and  quasi¬ 
static  loads.  The  usual  cartesian  tensor  indicial 
notation  and  summation  convention  are  used  in  some 
places  and  standard  scalar  notation  in  others. 

In  rectangular  cartesian  coordinates,  the  equili¬ 
brium  equations  are 

d  a,  +  f,  =  0.  (l) 

m  km  k 

The  body  force  f,  will  be  negligible  in  the  experiments. 

The  a  are  the  components  of  the  stress  tensor  which 
km 

is  symmetric.  The  symbol  will  mean  partial  deriva¬ 
tive  with  respect  to  the  rectangular  cartesian  coor¬ 
dinate  xm,  throughout  the  following  discussion. 

The  components  of  the  strain  tensor  are 

related  to  the  displacement  vector  ufc  by 

0  ~  \  (ch  u  +  B  uv) .  (2) 

km  2  v  k  m  m  k' 


In  order  to  write  the  stress-strain  relations 
it  is  convenient  to  introduce  the  deviatoric  stress  s 
and  deviatoric  strain  e^  defined  by 

s km  =  CTkm  ‘  1/3  Grr  6km'  (3) 
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(4) 


e,  =e,  -  l/3  e  6,  • 

km  km  rr  km 

where  6,  denotes  the  Kronecker  delta.  Then,  for  an 
km  9 

isotropic  material 

skm  "  Jo  2G(t‘T)  ®kn,(T)  dT>  (5) 

!/3  orr  ■  J*  K(t-r)  €rr(r)  dT.  (6) 

Where,  the  dot  indicates  the  derivative  with  respect 
to  the  argument  indicated. 

The  inverses  of  equations  (5)  and  (6)  are 

2ekm  '  Jo  J(t-T)  dkn,<T>  dT-  <7> 

err  =  1/3  B ( t- t )  arr(T )  dT.  (8) 

The  functions  G(t)  and  K(t),  or  J(t)  and  B(t) 
characterize  the  material  and  have  the  names  ; 

J(t)  =  shear  creep  compliance* 

B(t)  =  hulk  creep  compliance  • 

G(t)  =  shear  relaxation  modulus* 

K(t)  =  hulk  relaxation  modulus. 

The  material  will  usually  he  supposed  homogeneous 
so  that  all  functions  are  independent  of  the  space  co¬ 
ordinates  . 

The  functions  are  related  hy  Volterra  integral 
equations : 


k 


(9) 


G(  t-  T  )  J(  T  )  dT  =  t, 

*  o 

\l  K ( t- t )  B(t)  dT  =  t.  (10) 

<j  0 

Given  J(t),  G(t)  can  be  calculated  by  solving  (9)  by 
the  method  shown  in  section  2.9;  etc. 

The  uniaxial  stress  field  which  is  approximated 
by  the  tension  test  is  of  special  importance.  If 
a  /  0  while  the  other  components  of  stress  are  zero, 
the  deviatoric  stresses  are 

S11  =  CTU'  3 22  =-1/;3  Gll’  s 33  =  "1//3  CT11* 

(11) 

and  the  others  are  zero.  From  equation  (4),  (7),  and 

(8), 

e1]L  =  D(t-r)  ct11  dT  ,  (12) 

e22  =  e33  =  '  Jo  v(t"T)  D(t_T)  °U(t)  dT*  (13) 

Where , 

D(t)  =  1/3  J(t)  *  1/9  B(t),  (14) 

v(t)  ,  >  ilil  ~  2  -  .  (15) 

6  J(t)  +  2  B(t) 

The  pair  of  functions  D(t)  and  v(t)  also 
characterize  the  material  and  have  the  names: 
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D(t)  -  tensile  creep  compliance, 
v(t)  -"Poisson’s  Ratio”  for  creep. 

This  definition  for  v  leads  to  =  -v  ex  for  the  uni¬ 
axial  creep  test.  This  definition  is  not  a  unique 
generalization  of  the  idea  of  Poisson’s  ratio  for 
elastic  materials . 

Solving  equations  (14)  and  (15)  we  have 

J  =  2  (l+v)  D,  (16) 

B  =  3  (l-2v)  D.  (17) 

Thus,  the  bulk  creep  compliance  and  the  shear  creep 
compliance  and,  consequently,  the  relaxation  moduli 
can  all  be  calculated  from  the  functions  v(t)  and  D(t) 
occurring  in  the  equations  describing  the  tensile  test. 
The  inverse  of  (12)  is 

an  =  k  E(t-T)  dT*  (i8) 

The  new  function  is  related  to  the  tensile  creep  com¬ 
pliance  by  the  Yolterra  integral  equation, 

ft  D(t-T)  E(t)  dT  =  t,  (19) 

J  o 

and  has  the  name : 

E(t)  -  tensile  relaxation  modulus. 

The  complete  solution  to  a  problem  in  linear 
viscoelasticity  is  determined  by  equations  (l)  -  (4) 
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and  a  pair  of  the  stress-strain  relations  such  as 
(5)  and  (6). 

The  solution  must  satisfy  certain  boundary 
conditions:  Usually  the  displacement  vector  is  given 

over  part  of  the  boundary  while  the  stress  vector  Fk 
is  given  over  the  remainder.  The  stress  vector  is 
related  to  the  stress  tensor  by 


F,  =  ct,  n  , 

k  km  m 

where  nm  are  the  components 
surface . 


(20) 

of  the  unit  normal  to  the 


2.2.  Elasticity. 

The  theory  of  linear  elastic  materials  may  be 
viewed  as  a  special  case  in  the  theory  of  viscoelas¬ 
ticity  when 


J(t)  =  JQh(t) 

and 

B(t)  =  Boh(t)  , 


where  h(t)  is  a  unit  step  function: 


h 


t<  0 
t>  0  * 


(21) 

(22) 


(23) 


Then 

G  =  Gq  h(t), 
K  -  Ko  h(t) , 


(24) 
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(25) 


D(t)  =  Dq  h(t), 
E(t)  =  Eq  h(t). 


The  elastic  constants  G  ,  E  . 

o’  o’ 

v  ,  and  K  are 

0 9  0 

rel ated  : 

G°  2  (l*vc) 

(26) 

E 

K  =  - ^ -  . 

3  (l“2vo) 

(27) 

The  field  equations  become 

a  a.  +  f1  =  0  , 

m  km  k  9 

(28) 

ekm  2  ^k  Um  +  ^m  Uk^ 

(29) 

Skm  =  Gkm  "  Grr  6km' 

(30) 

e,  =  e,  “  l/3  e  6n  , 

km  km  '  rr  knr 

(31) 

s ,  =  2  G  e  , 

km  0  km 9 

(32) 

l/3  a  =  K  e 

rr  0  rr 

(33) 

Thus,  time  plays  the  role  of  a  parameter  and  does  not 
appear  explicitly  in  the  field  equations. 
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2.3®  Proportional  Loading; » 

An  important  class  of  problems  are  those  for  which 
the  stress  vector  on  the  surface  has  the  form 


Fk  =  Fk  f(t) 


(34) 


The  superposed  har  will  denote  a  function  of  the  space 
coordinates  only.  This  is  termed  proportional  loading. 
Let  us  determine  under  what  conditions  the  stress  field 


has  the  form 


°km  ‘  akm  f(t)> 


(35) 


where  a  can  be  determined  from  the  theory  of  elasti 
km 


city. 


The  boundary  conditions  (20)  give 


aVmn  =  F,  . 
km  m  k 


(36) 


In  the  absence  of  body  force^  the  equilibrium  equations 


(1)  give 


S  ct,  =  0  . 
m  km 


The  deviatoric  stresses  are 


(37) 


s,  =  s ,  f(t), 
km  km  ’ 


(38) 


where 


s km  =  CTkm  '  1/3  arr  6km  * 


(39) 


The  stress-strain  relations  (7)  and  (8)  give 

=  e,  g(t)  ,  (40) 


9 


where 


e  =  e  k(  t )  ; 
rr  rr  '  '  ’ 


=  1/ G  s.  , 

'  0  km  9 

(41) 

1/3K  a 

0  rr  , 

(42) 

G  J(t-T) 

f(T)  dT  , 

(43) 

rt  .  . 

K  B(t-r) 

U  t)  * 

f  (  T  )  dT  . 

(44) 

The  constants  G  and  K  are  arbitrary  constants, 
o  o 


e,  =  e .  +  1/3  6.  e 

km  km  km  rr 


Then 


1/2  (a,  u  +  a  u  ), 

'  x  k  m  m  k'  ’ 


(45) 

(46) 


Uk  =  uk  g(t)  »  (46) 

only  if  g(t)  =  k(t)  or  =  0;  that  is,  only  when 

B(t)  is  proportional  to  J(t)  by  a  constant  factor. 

From  equations  (16)  and  (17), 

B(t)  =  3/2  1~2  J(t)  .  (47) 

l+v(t) 

Thus  B(t)  is  proportional  to  J(t)  only  when  v  is  con¬ 
stant.  Further,  =  0  only  when  v  is  a  constant 

equal  to  1/2. 

If  v  is  constant,  then  the  stress  and  displace¬ 
ment  are  determined  by  (35)  and  (46)  where  the  barred 
quantities  satisfy  the  equations  of  elasticity,  pro¬ 
vided  the  boundary  conditions  of  the  viscoelastic 
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problem  agrees  with  (34) >  and  (46). 

In  fact  v  is  never  constant  for  any  material.  But 
it  frequently  happens  that  the  variation  in  v  with  time 
is  slight  or  the  effect  of  changes  in  v  on  the  solution 
to  the  elasticity  problem  is  slight.  In  these  cases 
the  stress  field  under  proportional  loading  can  be 
approximated  by  the  elasticity  solution. 


2.4.  Plane  Stress. 

If  the  stress  components  0^-3  =  0  >  "the  equations 
of  viscoelasticity  can  be  simplified  to  the  following: 


Vae  =  0  ’ 


ea3  1//2  (Ve  +  ’ 


s  «  =  a.. o  -  1/3  a....  6 


a.6 


a8 


rr  a8  ’ 


ea8  ea6  "  err  6a8  ’ 


2e 


a6 


rr 


J(t-r)  sa0(T)  dT  » 
1/3  B(t-T)  Cfaa(T)dT  • 


(48) 

(49) 

(50) 

(51) 

(52) 

(53) 


The  Greek  indices  have  the  range  1,  2  and  all  quantities 
are  functions  of  xa  and  t. 

The  equations  (48)  -  (53)  define  a  plane  stress 
problem.  There  is  an  important  class  of  plane  stress 
problems,  for  which  the  solution  is  independent  of  the 
material  properties,  which  we  now  consider.  In  this 
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case,  the  stress  can  be  calculated  from  the  corres¬ 
ponding  elasticity  solution. 

Equations  (4S)  are  satisfied  by  introducing  the 
Airy  stress  function  F: 

a  0  =6  &  a*  F  0  (54) 

a6  ay  p5  y  6 

Where  is  the  permutation  symbol  (£^  =  “  ^1  =  ■*** 

£11  =  ^22  =  * 

Eliminating  the  displacements  from  (49)  gives 


5  5ceo-0. 

°a6  8  t  t  6  a  8 


(55) 


This  is  the  compatibility  relation. 

Substituting  (54)  into  (50)  -  (53)  and  the  result 
into  (55)  leads  to 


v  F  =  3  a  F  =  0 . 

a  a  0  3 


(56) 


If  the  boundary  conditions  are  solely  on  stress,  the 
equations  (56)  and  the  stress  boundary  conditions  com¬ 
pletely  determine  the  stress  independent  of  the  material 
properties  and  the  stress  is  the  same  as  for  an  elastic 
material . 


2.5.  Short  Time  Solution. 

This  section  deals  with  a  suddently  applied  load 
and  the  solution  immediately  after  loading  (that  is, 
immediately  after  all  wave  propagation  effects  are 
completed ) . 


The  stress  and  strain  fields  are  related  hy  (5) 
and  (6).  For  a  step  change  at  time  zero, 

=  2G(0)  eta.  ’  (5?) 

1/3  arr  -  K(o)  €rr  .  (58) 

Thus,  constitutive  relations  are  identical  to  the  el¬ 
asticity  equations  with  elastic  moduli  equal  to  the 
values  of  the  relaxation  modulus  at  time  zero.  There¬ 
fore,  the  stress  distribution  at  the  instant  of  loading 
will  be  the  same  as  given  by  the  elasticity  solution 
for  shear  modulus  G  =  G (q)  an(l  bulk  modulus  K  =  K(o) 
or  Poisson  ratio  v  =  v(q)» 


2.6.  Long  Time  Solution. 

The  viscoelastic  materials  have  a  fading  memory. 
That  is,  the  stress  due  to  a  step  change  in  strain 
decreases  in  time.  The  result  is  that  when  loads  are 
maintained  at  a  constant  value,  an  equilibrium  state 
is  reached  such  that  no  further  change  in  stress  or 
strain  occurs . 


By  (5)  and  (6)^ 

S,  (<»)  =  G(co)  e,  (eo)  , 


(59) 

(60) 


km'"'  '  km 

1/3  CTrr(°°)  =  K(»)  5rr(°°) 

The  constitutive  equations  again  coincide  with  linear 
elasticity.  In  this  case  the  moduli  are  equal  to  the 


13 


values  of  the  relaxation  functions  at  time  infinity. 
Consequently,  if  the  surface  loads  remain  constant, 
the  stress  state  will  eventually  approach  that  of  the 
elasticity  solution  for  shear  modulus  G  =  G(°°)  and 
Poisson  ratio  v  =  v(°°). 

2.7.  Temperature  Effects . 

Each  of  the  material  property  functions  is  tem¬ 
perature  dependent.  That  is,  if  the  tensile  test  is 
conducted  at  different  temperatures,  the  tensile  creep 
compliance  is  a  different  function  of  time.  In  many 
instances  this  temperature  dependence  can  be  described 
by  an  empirical  relation  involving  the  so-called  "time 
temperature  equivalence".  The  material  is  then  termed 
thermo- rheologi c ally  s imple . 

Let  G(t)  be  the  relaxation  modulus  at  constant 
temperature  T  and  G(t)  be  the  relaxation  modulus  at 
some  reference  temperature  Tq : 


Ik 


Suppose  that  the  effect  of  temperature  on  the 
short  time  modulus  and  long  time  modulus  is  such  that: 

G(0)  =  G(0),  <61) 

G(<»)  =  G ( * )  .  (62) 

In  this  case,  there  always  exists  a  pair  of  times  for 
which  G  and  G  have  the  same  value 

G(t)  =  G(5). 

The  time  %  is  called  the  reduced  time.  It  depends  on  the 
temperature  T  of  the  test  and  the  time  t  at  which  the 
ordinate  G  is  measured. 

Empirical  evidence  suggest  that,  for  isothermal 
deformation  of  many  materials, 

5  =  a(T)  t.  (64) 

The  function  a(T)  of  temperature  is  known  as  a  tem¬ 
perature  shift  factor  (in  the  sense  of  shifting  the 
curves,  if  the  figure  were  plotted  on  log  time). 

Similar  equations  may  he  found  to  describe  the  iso¬ 
thermal  values  of  the  other  material  property  functions, 
hut  the  shift  factor  is  usually  different. 

Isothermal  experiments  cannot  determine  the  nature 
of  the  stress-strain  relationship  for  tests  in  which 
the  temperature  varies  with  time  and  no  experimental 
evidence  is  available.  A  possible  generalization  which 
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is  compatible  with  the  isothermal  observations  was 

* 

proposed  by  Morland  and  Lee.  Their  idea  is  that  rate 
of  change  of  the  modulus  is  determined  by  the  instan¬ 
taneous  temperature. 

If  G( t )  is  the  relaxation  modulus  measured  at 
varying  temperature  T(t),  then  (63)  holds.  Thus 

dG  dG  d_§ 

dt  "  d?  dt 

The  fundamental  hypothesis  is  that,  for  varying  tem¬ 
perature, 

||  =  a  ( T )  .  (66) 

This  agrees  with  (64)  in  the  isothermal  case.  If  tem¬ 
perature  is  a  function  of  time  for  a  given  particle, 
and  the  material  is  homogeneous, 

?(t)  =  f  a(T(t))  dt  .  (67) 

J  o 

This  generalizes  (64)  to  the  case  of  varying  tempera¬ 
ture  . 

Now  suppose  the  same  temperature  variation  occurs 
but  the  strain  is  applied  at  time  t.  Let  F  (t)  be  the 
resulting  relaxation  modulus.  Again  assuming  that  the 
initial  and  final  values  are  independent  of  tempera¬ 
ture,  there  exists  a  function  §  (t)  such  that 

FT(t)  =  G(?t)  ,  (68) 

*  Trans.  Soc.  Rheology,  _4,  233-263,  (I960). 
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and 


dt 


dt 


The  fundamental  hypothesis  is 

d£ 

— —  =  a(T)  . 

dt 


(69) 


Since 


Ft(t)  =  G(0)  , 

(70) 

1 

-R 

11 

O 

• 

(71) 

Thus  , 

§  =  P  a ( T ( t ) )  dt  =  ?(t)  -  5(t)  . 

t  J  o 

The  stress  due  to  continuously  changing  strain  is  then 


.  .  =  r  2  G(5(t)  -  §(t))  e .  . ( t )  dT.  (V2) 

1 J  <•'  o  d 


This  relation  is  a  hypothesis  not  yet  conformed  by 
experiment  and  it  can  be  expected  to  give,  at  best, 
only  approximate  representation  of  the  effects  of 
temperature . 


2.8.  Volterra  Integral  Equations. 

The  various  functions  of  time  which  describe  the 
physical  properties  of  a  viscoelastic  material  are 
related  by  integral  equations  of  the  form: 
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t 


(73) 


D(t-r)  E(t)  dT 


The  functions  D(t)  and  E(t)  have  the  general  form 


The  function  E(t)  decreases  monotonically  toward  a 
constant  value,  while  the  function  D(t)  increases  mono¬ 
tonically  to  a  constant  value.  The  limiting  values  are 
relat  ed : 


E(0) 


D(0) 


(74) 


and 


E(co)  = 


(75) 


D  ( 03 ) 

A  fundamental  problem  of  viscoelasticity  is  to  determine 
one  function  when  the  other  is  known. 

Equation  (73)  has  been  exhaustively  studied  and 
various  numerical  procedures  are  known  which  will,  in 
principle,  allow  the  solution  to  be  approximated  as 
closely  as  desired.  One  may,  for  example,  divide  the 
time  axis  into  finite  intervals  and  make  use  of  direct 
approximation  of  the  integral  by  some  numerical  means 
such  as  the  trapezoidal  rule. 
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In  principle,  numerical  integration  leads  to  the 
exact  answer  as  the  size  of  the  intervals  tends  to  zero. 
Unfortunately,  the  various  possible  numerical  proce¬ 
dures  that  may  he  used  to  approximate  the  integral  are 
very  sensitive  to  errors  in  the  values  of  the  functions 
and  such  errors  are  always  introduced  by  round-off 
of  the  numbers.  Moreover,  the  known  function  is  fre¬ 
quently  determined  by  experiment  and,  owing  to  experi¬ 
mental  error,  may  not  be  such  a  smooth  function  as 
the  exact  value.  Any  such  irregularities  tend  to  mag¬ 
nify  enormously  if  the  numerical  procedure  is  not 
chosen  carefully. 

Several  numerical  procedures  have  been  tried  in 

order  to  determine  E(t)  when  D(t)  was  known  from  ex- 

,  * 

periments.  The  procedure  of  Hopkins  and  Hamming  seems 
to  be  the  best. 


J.  Appl.  Phys . ,  28,  906-909,  (1957). 
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3.  TRANSMISSION  OF  LIGHT. 

3.1.  Plane  Waves, 

Consider  a  phenomenon  which  varies  with  distance 
z  and  time  t  according  to  the  law 

a  cos  ( kz  -  ojt  +  e  )  ,  (1 ) 

where  a,  k,  uu,  e  are  constants.  Such  a  phenomenon  is 
referred  to  as  a  wave.  The  factor  ,fan  is  the  amplitude 
and  the  argument  of  the  cosine  function  is  the  phase. 

When  considering  the  solution  to  a  differential 
equation  for  a  real  variable,  it  is  often  convenient 
to  consider  the  variable  to  be  a  complex  number.  The 
real  and  imaginary  parts  of  the  complex  variable  are 
then  each  solutions  of  the  differential  equation  for 
the  real  variable. 

The  function  (1)  is  the  real  part  of  the  complex 
variable 

ae1CD  ,  (2) 

where  a  is  a  complex  number  and 

cp  =  kz  -  art  .  (3) 

Any  phenomenon  represented  by  (2)  is  also  called  a  wave. 

This  representation  of  a  wave  by  complex  varia¬ 
bles  can  be  extended  to  vector  valued  functions:  A 
complex  vector  will  mean  an  entity  for  which  the  real 
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and  imaginary  parts  are  real 
denoted  by  a  subposed  tilde, 
has  the  form 


vectors.  Vectors  will  be 
A  vector  valued  wave  then 

(4) 


cd  =  k  n  •  r  -  uut. 

The  complex  constant  vector  a_  will  be  called  the  ampli 
tude  vector.  The  real  unit  vector  n  defines  the  di¬ 
rection  of  propagation.  The  real  vector  t  is  the  po¬ 
sition  vector.  The  frequency  of  the  wave  is  tu.  The 
wave  number  is  k.  The  wave  speed  is  defined  as 


The  wave  length  is  defined  as 


Let  the  real  part  of  a  complex  quantity  be  de¬ 
noted  by  a  superscript  +,  the  imaginary  part  by  a  super 
script  -,  and  the  conjugate  by  a  superscript  *.  The 
following  terminology  is  used. 


a  =  0 

circularly  polarized. 

a  =  0  - 

linearly  polarized  , 

n  =  0  - 

transverse  . 

r>->  r^J 

If  a+  a”,  and  n  form  a  right-handed  system,  the  wave 
is  called  right-handed.  Only  right-handed  transverse 
waves  are  considered  in  the  following. 
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The  real  part  of  the  complex  wave  is  given  by 
R(E)  =  a+  cos  cp  -  a  sin  cp 


Let  g.  be  the  base  vectors  and 
'■'*'1 


a-  g..  ,  a 
1  ~ 


a .  g . 

1  ~i 


Then 


R(E)  =  Ql  cos  (cp  +  e1)  J|1 
+  Q2  cos  (cp  +  e2)  £2 
+  Q  cos  (cp  +  e^)  } 

where 


(5) 


and 


a . 

tan  e.  =  — -  (no  sum) 

a . 

i 

Thus  the  vector  valued  wave  may  be  regarded  as  three 
waves;  each  has  the  same  frequency  and  wave  number  but 
different  amplitudes  and  different  phases*, 

The  case  of  a  transverse  wave  propagating  along 
the  z-axis  will  be  of  particular  interest*,  Then  a  + 
and  a  lie  in  the  x-y  plane  and  Q0  =  0.  The  vector 
expression  (5)  represents  two  waves  of  different  am¬ 
plitude  and  phase. 

For  the  linearly  polarized  wave,  a  x  a  =  o •  This 
implies  that  a.  x  a  =0.  Thus,  a.  and  a  are  parallel*, 
In  that  case,  the  phase  is  the  same  but  the 
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amplitude  may  be  different: 

R(E)  =  Q1  cos (co  +  e)  +  Q2  cos(cp  +  e)  £2.(6) 

For  a  circularly-polarized  transverse  wave  pro¬ 
pagated  along  the  z-axis,  the  requirement  a.  •  a.  =  o 
implies  the  a+*  a  are  orthogonal  and  have  the  same 
length.  In  this  case,  the  amplitudes  are  the  same  and 
the  phase  of  the  two  component  waves  differ  by  tt/2: 

R(E)  =  Q  cos (cd  +  e)  g1  +  Q  sin(cp  +  e  -  %)  g2#(7) 
3.2.  Dielectrics. 

This  section  is  devoted  to  the  theory  of  bire¬ 
fringence  based  upon  the  description  of  light  as  an  el¬ 
ectromagnetic  wave  propagating  in  an  anisotropic  di¬ 
electric.  The  applications  will  be  limited  to  quasi¬ 
static  deformations  so  that,  during  the  time  of  passage 
of  an  electromagnetic  wave,  the  displacements  of  the 
material  are  negligible.  The  equations  of  a  dielectric 
at  rest  can  therefore  be  used.  Furthermore,  since  only 
the  relative  intensity  of  monochromatic  light  will  be 
observed,  the  absorption  and  dispersion  will  be  neglected. 

The  equations  of  an  anisotropic  dielectric  then 
have  the  following  form: 
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curl 

S  “  Tt  =  0  • 

(8) 

curl 

bB 

I  =  0  • 

(9) 

D.  = 

e  K.  .  E  .  . 

(10) 

l 

o  1  J  J 

H  = 

—  B  . 

(11) 

u  ~ 

*0 

The  fields  D  and  B  are  also  subject  to  the  relations 
div  B  =  0  , 
div  D  =  0  ; 

but  these  will  be  satisfied  trivially  by  the  solution 

which  will  be  obtained  for  (8)  -  (11). 

The  tensor  K.  .  will  be  called  the  dielectric 
1 3 

tensor.  It  is  regarded  in  this  chapter  as  known 
and  constant.  The  term  "anisotropic " ,  as  used  here, 
merely  means  that  the  dielectric  tensor  does  not 
reduce  to  a  scalar  multiple  of  the  unit  matrix.  The 
and  are  fundamental  electromagnetic  constants. 

At  a  surface  with  unit  normal  vector  v,  the  fields 
are  subject  to  the  jump  conditions: 


O 

II 

i — i 

i _ i 

X 

>1 

X 

1 — I 

1 _ 1 

1! 

o 

(12) 

It 

• 

1 — 1 

1 _ 1 

II 

o 

* 

V.  [B]  =  0  . 

/-S-/ 

(13) 

It  will  turn  out  for  the  problems  considered  that  the 
relations  (13)  are  identically  satisfied  when  the 
relations  (12)  are  satisfied  and  consequently  they  are 
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not  used  here. 


The  index  of  refraction  N  is  defined  hy 


W  9 


(19) 


where 


(20) 


u  e 

o 


Eliminating  b  and  d  gives  the  equations 


a.  =  -  Li  c  N  K.  .  e  .  n  h  . 
l  ij  jrs  r  s  ’ 


(21) 


(6^  -  N2  6^'rm  K"1  n  n  )  h.  =  0 
l  ipq  qm  p  r  j 


(22) 


The  relation  (22)  represents  three  homogeneous 
equations  for  the  components  of  the  real  or  the 
imaginary  parts  of  h.  Thus,  a  wave  exists  only  when 
the  determinant  of  the  coefficients  is  zero: 


det  =  l  -  6S^  n  n,  K-,  ^  + 

ah  a  t  hs 


(23) 


4  !  /cst  *ik  ,ik  cStx  -1-1 

r-5*  (6,  6  -6,  6  )n  n,  n.  n  3L  K  .  =  o 

1  2  v  ha  me  ha  me/  a  t  k  e  hs  mi 


The  6-  symbols  are  the  generalized  Kronicker  delta. 

Since  N  >  o,  there  are  exactly  two  indices  of 
refraction,  and  N^,  determined  hy  the  roots  to  (23). 
If  Kj  .  is  symmetric  and  /  N^,  the  two  corresponding 
solutions  of  (22)  are  orthogonal  vectors.  Consider  the 
coordinate  system  with  the  z  axis  along  n  and  the  x,  y 
axes  along  the  directions  determined  by  the  solutions 


2  6 
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The  ,  K2,  K3  are  the  principal  values  of  the  tensor 

K.  ..  Similar  results  hold  when  n  is  normal  to  the  first 
1  J 

principal  direction. 

In  every  case  the  real  and  imaginary  parts  of 
each  field  are  determined  by  the  same  relations.  They 
are  therefore  parallel  and  each  field  is  a  linearly 
polarized  wave  in  the  sense  of  section  3.1.  The  field 
H  is  also  a  transverse  wave. 

The  direction  of  a  is  called  the  direction  of 
vibration.  The  direction  of  h  is  called  the  direction 
of  polarization.  It  will  be  most  convenient  in  the 
following  sections  to  select  the  vector  h  as  fundamental. 
The  vector  a  is  then  determined  by  one  of  the  relations 
(24)  -  (27). 


3.3.  Normal  Incidence. 

The  solution  of  equations  (8)  to  (11)  subject  to 
appropriate  boundary  conditions  gives  the  laws  of  the 
reflection  and  refraction  of  light.  These  results  are 
usually  given  in  the  form  of  geometric  constructions 
which  are  very  tedious  to  follow  or  else  in  approximate 
forms  for  which  the  nature  of  the  approximation  is  not 
clear.  The  following  sections  give  an  analytical 
expression  to  the  optical  laws  which  will  be  needed 
in  photomechanics. 
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Consider  a  slice  of  the  dielectric  which  is 
assumed  to  he  homogeneous,  i.e.,  K.  .  is  constant. 
Choose  the  coordinate  axes  so  that  the  surfaces  of 
the  slice  are  z  =  o  and  z  =  h. 

If  a  plane,  transverse  wave  is  incident  normal 
to  one  of  the  surfaces,  part  of  the  wave  is  reflected 
and  part  is  transmitted.  A  similar  reflection  and 
transmission  occurs  at  the  second  surface.  There  are, 
therefore,  five  electromagnetic  waves  having  the 
following  directions : 


A  solution  to  the  field  equations  is  required 
when  the  incident  wave  is  a  plane  transverse  wave 

Si  =  Si  elCPl>  cpi=ko2i’  £  "  UJt* 

If  h^  is  parallel-  to  one  of  the  secondary  prin¬ 
cipal  directions  of  K. .  in  the  plane  z  =  constant,  the 
fields 

H  =  h  eicp2  m  =  k  n  .  r  -  tut  , 

7  2  o  ~  7 

H0=  h„  elcc)3  cp  =  k  n  •  r  -  art  , 

^  J>  ^  j>  7  ~  7 
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art  , 
out  , 


have  h  ,  h_,  h,,  and  h,.  parallel  to  that  principal 
direction. 

There  are  two  fundamental  cases,  according  to 
whether  h  is  parallel  to  one  or  the  other  principal 
directions.  The  internal  fields  in  each  case  are  de¬ 
termined  by  (24 )  or  (25).  The  solution  for  an  ar¬ 
bitrary  wave  is  a  linear  combination  of  these  two 


cases. 

In  each  case,  the  jump  conditions  (12),  applied 
at  each  face,  lead  to  the  solution 


(28) 

(29)  . 


where  2N 


and  N  is  given  by  (24)  or  (25). 


(30) 

(3D 
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3.4.  Oblique  Incidence. 

Consider  the  case  when  the  normal  to  the  surface 

of  the  slice  of  dielectric  is  a  principal  axes  of  K.  .. 

1  J 

Choose  the  coordinates  along  the  principal  axes  of  K... 

1  j 

The  instance  when  the  wave  propagated  in  the  dielectric 
is  perpendicular  to  the  y-axis  will  now  be  considered. 

The  secondary  principal  axes  of  then  coincide 

with  the  y-axis  and  the  normal  to  the  y-axis.  The  two 
waves  given  by  (26)  and  (27)  are  then  possible. 

Consider  first  the  case  when  the  incident  wave 
has  h  parallel  to  the  y-axis.  The  reflected  and  trans- 
mitted  waves  turn  out  to  all  have  h  parallel  to  the 
y-axis.  The  five  fields  have  the  directions  shown: 


Applying  the  jump  conditions  (12)  to  these  waves 
and  using  the  solution  (26)  leads  to  the  relation 
between  incident  and  transmitted  waves : 

Hc  =  C  ei(6  "  koh  COS  tt)  H  °2) 
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and 


(33) 

(34) 

(35) 


y  =  kh  cos  8 

N  cos  6 
K^cos  a  * 


(36) 

(37) 


and  Kx,  K2,  K^,  are  the  principal  values  of  K^.. 

Consider  next  the  case  when  the  incident  wave  has 
h  normal  to  the  y-axis.  Then  the  solution  (27)  applies. 
Using  the  jump  conditions,  the  relation  between  incident 
and  transmitted  waves  is  again  found  to  have  the  form 
given  above  but  with  f  and  N  given  by  new  formulae: 


N  cos  8 
cos  a 


(40) 


N  =  N2  =  /K2 


(41) 
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Every  incident  wave  normal  to  the  y-axis  can  be 
represented  as  a  linear  combination  of  these  two  cases. 

Similar  results  hold  when  the  incident  wave  is 
normal  to  the  x-axis. 

When  the  wave  is  along  the  z-axis,  the  results 
coincide  with  equation  (30)  and  (31). 


3o5.  Slightly  anisotropic,  transparent  materials. 

The  materials  used  in  photomechanics  are  trans¬ 
parent  isotropic  materials  for  which  straining  causes 
a  small  change  in  dielectric  properties.  Let  Nq  denote 
the  index  of  refraction  in  the  undisturbed  state  and 
consider  the  propagation  of  light  through  a  slice  of 
thickness  h.  The  dielectric  tensor  can  be  written 


K.  .  =  N2  (6.  .  +  ^-r  P.  .)  , 
V  n  11  "'V'  1  1  7 


1  J 


1  J 


ujh 


(42) 


where  P.  .  depends  on  the  strain  history  and  is  zero 
i  J 

in  the  undisturbed  state. 

The  exact  relations  of  sections  3.3.  and  3.4.  for 
propagation  of  light  through  a  slice  will  now  be  in¬ 
vestigated  for  the  case  when 

a)  the  material  is  slightly  anisotropic,  i.e., 
when  P.  .  is  of  order  unity,  and 

b)  the  material  is  nearly  transparent,  i.e., 
when  Nq  is  of  order  unity. 


There  are  two  kinds  of  approximation  involved. 

The  first  is  that  the  second  term  in  (1)  is  small 
becaus e 

€  =  —  «  1  .  (43) 

xh 

The  parameter  e  is  a  non-dimensional  characteristic 
number  of  the  problem,  usually  having  a  magnitude  about 
10-5.  Therefore  terms  of  higher  order  in  e  can  be 
neglected.  The  second  is  that  the  components  K^.  are 
of  order  unity.  In  fact  Nq  is  about  1.5,  but  the  error 
involved  by  the  second  approximation  does  not  exceed 
a  few  per  cent,  especially  when  only  relative  retarda¬ 
tion  is  observed. 

Consider  first  the  results  for  normal  incidence. 
The  index  of  refraction  is  given  by  (24)  or  (25). 

Thus  , 

N1  =  J  llj  (1  *  t  P1;L)  =  N0(l  *  J  :  pu)>  U4) 

“2  =  J  "o  (1  *  £  P22>  :  No(1  *  *  £  P22>-  U5) 

Equations  (44)  and  (45)  are  of  the  form 

N  =  N  +  e6  .  (46) 

o 

Thus 

kh=— =N  — +6.  (4?) 

e  o  e 
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For  simplicity,  suppose  that  e  is  chosen  such 


that 

the  large  number  N  — 

O  G 

is 

an  integer  multiple 

of 

2n. 

Then 

sin  kh  =  sin 

(N  -  + 

v  OG 

6) 

=  sin  6  , 

(48) 

cos  kh  =  cos 

(N  -  + 
v  oe 

6) 

=  cos  6  , 

(49) 

The 

relation  (31) 

becomes 

,  _  *  1+N§ 

tan  A  =  - 

tan  6 

(50) 

2N 

o 

for  slightly  hiref ringent  materials. 

2 

If  N  is  1.5,  1  +  Nn  =  1.08  . 

°  2N 

o 

As  approaches  unity, this  factor  approaches  unity. 
Thus,  in  the  two  cases  (44)  and  (45), 


91  =  ±  No  P11 


(51) 


and 

P22  '  <52> 

for  nearly  transparent  materials. 

The  relation  (30)  "becomes 
2N 

o 

C  =  /  p  p  n  p  9  (53) 

J  4  N  cos  p  +  (1+N  )sin  0 
V  O  o 

for  slightly  hiref ringent  materials. 
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The  absolute  value  of  C  has  the  bounds 


<  C  <1. 


The  value  is  therefore  nearly  unity  for  nearly  trans¬ 


parent  materials; 


C  =  1  . 


(54) 


Secondly,  consider  the  problem  of  oblique  inci¬ 
dence  solved  in  section  3.4.  From  (42), 


K1  =  ( !  +  e  ?±)  . 


(55) 


Thus , 


K11  =  No2  (1  "  e  Pl}  ‘ 


Similar  expressions  hold  for  K 2  and  . 

The  relations  (3?)  to  (39)  lead  to 


sin  8  =  -jj-  sin  a. 


N  cosa 
o 


Nx  =  No  (1  +  i  P1  +  i  P3  tan  p) 


(56) 


(57) 


Thus  , 

N  cos  P 

Y  =  - " -  +  6  , 
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where 


6  =  \  Nq  cos  0  (P1  +  tan  0)  . 

Choosing  e  so  that  (Nq  cos3)/e  is  a  multiple  of  2tt, 
(35)  becomes 


(58) 


1  +  f 

tan  (1  =  — —  ■■  tan  6  . 


(59) 


If  the  angle  of  incidence  varies,  for  example, 
from  zero  to  45°  and  Nq  =  1.5, 


0.667  <  f  <  0.833  . 


In  this  case 


1.01  <  -  <  1.08. 


Therefore,  only  a  small  error  is  made,  if  (59)  is  re 
placed  by 


0  =  6 


(60) 


and  (34)  is  replaced  hy 


C  =  1  . 


(61) 


If  equations  (40)  and  (41)  are  used,  then  (60) 


and  (61)  again  hold  "but  with 


2  cos  B 
in  place  of  ( 58) . 


(62) 
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The  two  waves  transmitted  at  oblique  incidence 


therefore  have,  approximately. 


o 

-1 


0 

-2 


cos 


cos  3 


6 


} 


"o  SiA 

2  cos  3  } 


i 


(63) 

(64) 


when  the  material  is  slightly  anisotropic  and  nearly 
transparent 0  These  equations  reduce  to  (51)  and  (52) 
for  normal  incidence. 


3.6.  Approximate  Optical  Relations. 

The  exact  expressions  of  sections  3.3  and  3.4 
for  transmission  of  light  through  a  slice  of  dielectric 
are  rather  complicated  because  of  the  losses  due  to 
reflection  at  the  entering  and  exiting  surfaces.  It 
was  shown  in  section  3*5  how  the  formulae  for  light 
transmitted  through  a  slice  could  he  approximated  when 
the  material  is  nearly  transparent  and  only  slightly 
anisotropic.  The  approximate  relations  are  now 
summarized . 

Consider  the  case  of  normal  incidence  when  the 
normal  to  the  surfaces  of  the  slice  is  a  principal 
axis  of  the  dielectric  tensor.  Let  £  he  a  unit  vector 
along  the  first  principal  axis  (with  index  of  refraction 

Nn  and  wave  number  )  and  q  he  a  unit  vector  along 

1  ±  ~ 


59 


the  second  principal  axis  (with  index  of  refraction 
and  wave  number  k^)  of  the  dielectric  tensor.  Let  the 
incident  wave  be 


IT  ,  ICO 

H  T  =  h  e  , 


(65) 


hT  =  hn  P  +  h„  q. 


(66) 


where  h^  and  h^  may  be  complex  constants.  The  trans¬ 
mitted  wave  is  given  by  (28)  applied  to  each  component 
Using  the  approximate  expressions  (51),  (52)  and  (54), 
the  transmitted  wave  is 


H„  =  h  elcp  , 
~jL 


,  ,  i(0n  -  k  h)  ,  i (0  -k  h) 

hm  =  h_  p  e  v  1  o  +  h.  q  2  o  , 

~±  1  ~  2  ~ 


ei  =  i  No  Pi  > 


02  =  ^  No  P2  * 


(67) 

(68) 

(69) 

(70) 


In  the  case  of  oblique  incidence  at  angle  a  but 
normal  to  the  second  principal  axis  of  the  dielectric 
tensor,  the  transmitted  wave  given  by 


ei(Bl-koh  oosa)*  hz2  ei(e2'koh  003  a),(71) 


N0  2  No  sin^p 

ei  =  t  pi  cos  8  +  t  p3  Tzn 
e  -l2_ 

d2  2  cos  0  ' 

1 

sm  0  =  —  sm  a  . 


(72) 

(73) 

(74) 


40 


The  exact  expressions  for  transmission  of  light 
may  also  be  approximated  by  simple  relations  in  the 
case  of  certain  special  materials.  Two  examples  - 
the  polarizer  and  the  quarter-wave  plate  -  are  of 
particular  interest. 

An  ideal  quarter-wave  plate  is  a  slice  of  di¬ 
electric  for  which  waves  with  h  along  the  first  prin- 

/■v/ 

cipal  axis  are  not  retarded  while  waves  with  h  along 
the  second  principal  axis  are  retarded  by  ninety  degrees. 
Thus  (68)  becomes 


=  h,  P 

i  ~ 


h2  e 


“  1 


rr/  2 


(75) 


The  real  part  of  the  incident  wave,  according  to  (5) 
is 

R(Hj )  =  Q1  cos(cp  +  £-]_)£  +  Q2  cos  (<p  +  e  2)  q.(76) 


The  real  part  of  the  transmitted  wave  is 

R(Ht)  =  Qx  cos(cp  +  Gq)  P  +  Q2  cos(cp  +  e2-n/2)  q 

(77) 


A  situation  of  particular  interest  occurs  when 
the  incident  wave  is  a  plane  polarized  wave  with  h 
making  an  angle  of  45°  from  jd  toward  Then  Q-^  =  Q2  = 

Q  and  eq  =  e2  =  e#  The  'transnii‘t'ted  wave  is,  according 
to  (7),  circularly  polarized. 


An  ideal  polarizer  is  a  dielectric  for  which 
the  first  principal  index  of  refraction  is  nearly  unity 
and  the  second  is  very  large.  From  (30),  it  is  seen 


that  when  h  is  parallel  to  the  first  principal  axis 

/-'V 

C  =  1  and  when  h  is  parallel  to  the  second  principal 
axis  C  =  0.  Given  an  incident  wave 


H 


(78) 


the 


transmitted  wave  is 
Ht  =  (h  •  £)  £  elcp. 


(79) 


The  real  part  of  the  incident  wave,  referred  to 
arbitrary  x-y  axis  in  the  plane  of  the  slice,  is 


R(H];)  =  Q1  cos  (cp  +  e]L)  +  Q2  cos(cp+e2)  g2- 


The  real  part  of  the  transmitted  wave  is 
R(Ht)  =  Qx  cos  a  cos(cp  +  e1)  £ 

+  Q2  sin  a  cos  (cp  +  e^)  £ 


(80) 


when  a  is  the  angle  between  the  axis  of  the  polarizer 
and  the  x  “  axis . 


3.7.  Plane- Polaris  cope. 

The  term  polariscope  is  a  general  expression  re¬ 
ferring  to  any  of  various  possible  arrangements  of 
polarizer  and  quarter-wave  plates  in  order  to  view  a 


nearly  transparent  dielectric.  In  the  plane-polariscope, 
light  is  passed  through  a  polarizer,  and  then  through 
a  second  polarizer  (called  the  analyzer).  The  trans¬ 
mitted  light  exhibits  interference  fringes  which  are 
now  investigated. 

Let  the  unit  vector  £  be  directed  along  the  axis 
of  the  polarizer.  Then,  by  (79),  the  light  from  a 
monochromatic  source  which  is  transmitted  by  the  pola¬ 
rizer  is  such  that 

H  =  A  £  elcp.  (81) 

The  origin  of  coordinates  can  be  chosen  so  that  A  is 
a  real  number. 

Choose  the  z-axis  normal  to  the  dielectric  and 
let  the  x-y  axes  be  directed  along  the  secondary  prin¬ 
cipal  axes  of  the  dielectric  tensor.  Only  the  case  of 
normal  incidence  will  be  considered,  but  identical 
results  hold  for  the  oblique  incidence. 

Let  C  be  the  angle  between  £  and  the  x-axis: 

x 

The  wave  incident  on  the  dielectric  is  such  that 
H  =  (a  cos  C  gp)+(A  sin  C  £2)  elcp. 
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(82) 


By  (68),  the  transmitted  light  is  given  by 

n  f  a  r  i  ( 0,  -  k  h ) 

H  =  |A  cos  C  e  ^  °  M, 

,  ■  ,  -?(&,  -  k  h)  1  icp 

+  A  sm  Q  e  4  o  g2  j  e  .  (830 

The  light  is  next  passed  through  the  analyzer. 

Two  cases  are  of  special  interest. 

Case  1:  Crossed  Polari zer-Analyzer . 

Let  the  axis  of  the  analyzer  be  described  by  the 
unit  vector  q.  Consider  the  case  when  q  makes  an  angle 

r-O  '■'•v 

Q  +  rr/2  with  the  x-axis: 


The  light  transmitted  by  the  analyzer  is  then,  by  (79), 


H 


f  .  -  .  -  i  (0,  -  k  h) 

=  j-  Acos  Q  sm  C,  e  v  1  o 


+  A  sin  Q  cos  Q  e 
H  q  . 


i  (0,  -  %  e 


icp 


(84) 


The  real  part  is 

R(H)  =  A  sin  2Q  sin  rnr  cos  (co  -  kQh  -  e), 

where 


2tt 


(85) 


(86) 
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The  variable  n  is 


and  e  is  determined  by  0^  and  0^. 
called  the  fringe  order. 

It  can  be  seen  from  (85)  that  there  will  be  zero 
light  intensity  at  those  points  for  which  the  principal 
axes  of  the  dielectric  make  an  angle  of  zero  or  90° 
with  the  axis  of  the  polarizer.  The  locus  of  all  such 
points  form  lines  on  the  model  called  isoclinic  lines. 

There  will  also  be  zero  light  intensity  where  the 
fringe  order  n  is  an  integer.  The  locus  of  all  such 
points  form  lines  on  the  model  which  are  called  iso- 


chromatic  lines. 

If  the  loss  of  light  by  reflection  is  not  neglec¬ 
ted,  it  is  found  that  the  isochromatics  are  only 
approximately  given  by  integer  fringe  orders  and  the 
intensity  is  not  zero.  But,  this  approximation  seems 
to  be  adequate  for  nearly  transparent  materials. 

Case  2:  Parallel  Polarizer- Analyzer . 

Let  the  axis  of  the  analyzer  have  the  same  di¬ 
rection  as  the  polarizer,  i.e.  £  =  £  .  Then,  by  (79), 
the  light  transmitted  by  the  analyzer  is 


H  = 


f,  2r  1(0,  -k  h) 

|  A  cos  C  e  1  ° 

.2-  i(ea  -k  h)}  elcp  p 
+  A  sm  Q  e  v  2  o')  * 


(87) 


H  £,  • 
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The  real  part  is  given: 


R(H)  =  A  Q  cos  (cp  -  kQh  +  e)  ,  (88) 

where 

02  =  cos +  sin^  +  2  cos2^  sin2£  cos  2nirt  (89) 

The  minimum  values  of  Q  occur  when  n  is  an  odd 
multiple  of  The  isochromatic  lines,  in  this  case 

connect  all  points  for  which  n  is  \9  3/2,  etc.  They 
are  referred  to  as  half-order  fringes. 


3*8.  Circular- Polaris  cope . 

The  isoclinic  lines  and  the  isochromatic  lines  are 
superposed  in  the  case  of  the  plane  polariscope.  In 
order  to  observe  the  isochromatic  pattern  more  easily, 
a  circular-polar is  cope  is  used.  The  circular-polaris cope 
consists  of  polarizer,  quarter-wave  plate,  model,  quarter 
wave  plate,  and  analyzer.  As  will  now  he  shown,  the 
isoclinic  lines  do  not  appear  on  the  model. 

Let  axis  of  the  polarizer.  Let  and 

q  he  the  axis  of  the  quarter-wave  plate.  Only  the  case 
when  2i  m&kes  an  angle  of  45°  between  and  ^  will  he 
cons idered : 
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The  light  transmitted  by  the  polarizer  and  quarter-wave 
plate  combination  is,  according  to  section  3.6,  cir¬ 


cularly  polarized: 


£  "  A  (£?"  1  6  A  ^~l" 


i(<P+C) 


(90) 


Again  A  can  be  chosen  to  be  a  real  number. 

The  light  is  then  passed  through  a  nearly  trans¬ 
parent  dielectric.  Let  the  x-y  axis  be  principal 
axis  of  the  dielectric.  The  transmitted  wave  is  given 


by  (68): 


tt  a  i(cp  +  h) 
H  =  A  e  '  o 


S,  e  -  i  &■ 


(91) 


G-koh)  {(e1 


Different  possibilities  now  arise  depending  on  the 
orientation  of  the  second  quarter-wave  plate  and  the 
analyzer.  Only  four  will  be  considered. 

Case  1:  Quarter-Wave  Plates  Augment  and  the  Polarizer- 
Analyzer  are  Parallel. 

Let  and  give  the  axes  of  the  second  quarter- 
wave  plate  and  the  axis  of  the  analyzer.  Only  the 
case  when  makes  an  angle  of  45°  between  and  3>3 
will  be  considered.  Now,,  let  the  quarter-wave  plates 
augment  and  the  polarizer-analyzer  be  parallel. 
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i.e 


£3  =  £2'  %  =  £2’  £4  =  £l  : 


The  wave  transmitted  by  the  quarter-wave  plate  is 


.  i ( cp  +  C“  k  h)f  /  i 0 1  _  _i9 

H  =  A  e  v  b  o  M (e  1 


cos  C  "i  sin  Q)^ 


(-e^A  cos  Q  +  i  e^f  sin  Q 


(92) 


The  wave  transmitted  by  the  polarizer  is  given  by 


£  "  H  £4  ’ 

H  =  i  r  ei(cp*tC-V>>{ei0-  -e16*} 

The  real  part  is 

R(H)  =  Q  cos  (cp  +  2 C  -  kQh  +  e) 


(93) 

(94) 


(95) 


where 

Q  =  A  /?  sin  2mr  .  (96) 

Thus,  isochromat ic  lines  give  integer  fringes. 
The  isoclinics  do  not  occur. 
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Case  2: 


Quarter-Wave  Plates  Opposed.  Polarizer- 


Analyzer  Crossed. 

Consider  the 


cas  e 


when  £3  =  £2,  £3  =  -  £2: 


The  wave  transmitted  by  the  quarter-wave  plate  is 

H  =  A  ei(co  +  £  “  ^o^^Ci  cos  Q  e^®1  +  sin  Q  e5'®*  )q3  (97) 

+  (-sin  C  "  i  cos  Q  ei®a  )  p^j.. 

The  wave  transmitted  by  the  polarizer  is 
H  =  A  ^  ei(cD  +  2C  “  koh  +f  ){ei9^-eie^  }.  (98) 

The  real  part  is 

R(H)  =  A  /2  sin  2nrr  cos  (cd  +  2£  -  kQh  +  cj  +  e).  (99) 

In  this  case  also,  the  isochromatic  fringes  are  of 
integral  fringe  order. 

Case  3:  Quarter-Wave  Plate  Aligned.  Polari zer- Analyzer 
Crossed . 

The  first  two  cases  give  a  dark  field  in  the 
absence  of  a  model  and  yield  integral  order  fringes. 


By  effectively  rotating  the  analyzer,  a  light  field  is 
obtained  which  yields  half-order  fringes. 


Consider  the  case  =  £2*  The  wave 

transmitted  by  the  second  quarter-wave  plate  is  then 

H  =  =  A  +  ^  ^oh^|(cos  Q  e^^  -i  sin  Q  e^2  ) (100) 

+  (-i  sin  Q  eif^  +  cos  Q  ei®^  )q  } 

The  wave  transmitted  by  the  analyzer  is 

tt  1  i  (cp-k  h)  /  i0f  ^  ifb  v  *  /n  nn  \ 

H  =  A  -j2  e  o  '  (e  1  +  e  2  ).  (101) 

The  real  part  is 


R(H)  =  A  /2  cos  nn  cos  (cp-  kQh  +  e)Q  (102) 

Thus ,  dark  lines  correspond  to  prints  for  which  n  is 
an  odd  multiple  of  \ . 

Case  4:  Quarter-Wave  Plates  Oppose.  Polarizer- 
Analyzer  Parallel . 

To  = 

transmitted  by  the  analyzer  has  the  real  part: 

R(H)  =  A  /2  cos  nn  cos  (cp-  k  h  +  +  e)  (103) 

O  o 

The  isochromatic  lines  are  again  half-order  fringes. 


uonsiaer  tne  case 


£3 


£2 


me  wave 
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In  each  of  the  four  cases  the  phase  of  the  wave 


contains  a 
but  is  not 


factor  e  which  is  determined  by  0]L  and  @2 
listed  because  it  does  not  effect  the  fringe 


pattern 


4 


BIREFRI NGENCE. 


4.1.  Rheo-Optic  Constitutive  Relations. 

The  materials  used  in  photo-mechanics  are  such 
that  the  dielectric  tensor  is  a  function  of  the  history 
of  strain.  The  material  is  isotropic,  in  the  sense 
that  response  of  a  sample  cut  from  a  sheet  is  inde¬ 
pendent  of  the  orientation  of  that  sample.  For  small 
displacements,  the  relation  will  be  linear: 


N^6.  .  +  6.  .  I**  A(t-T)e,,  (t) 

o  ij  ij  J0  '  ' 

+  J*J  2B ( t-  t  )  e i  j  (  t  )  dT 
In  terms  of  the  mean  and  d 


dT 


eviatori 


c 


parts  , 


(i) 


K.  .-  1/3 
1 J 


IJ 


K 


rr 


3 


the  relation  (1)  has  the  form: 

Krr  =  3  No+  Jo  C3A( t- T ) +  2B(t-T) ]err (t)  dr,  (2) 

Rij  =  Jo  2B ( t- t ) j ( t )  dT  .  (3) 

Using  the  stress-strain  law  of  chapter  1,  the  entirely 
equivalent  expression  for  the  dielectric  tensor  in 
terms  of  stress  history  is  found: 
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(4) 


Krr  =  3No  +  Jo  0(t-T)arr(T)  dr, 

R.  .  =  f*  D(t-r)s.  .(T)  dT,  (5) 

1  j  J  o  1  j 

where 

C(t)  =  J*J[3A(t-T)  +  2B(t-T)]  B(t)  dT  ,  (6) 

D(t)  =  B(t-T)  J(t)  dT  .  (7) 

The  deviatoric  part  is  especially  important  for 
the  interpretation  of  the  fringe  order.  The  inverse 
of  (5)  will  be  important.  It  is 

sij  =  Jo  E(t-T)  LjCt)  dT  ,  (8) 

where  E(t)  and  D(t)  are  related  by  the  integral  equation 
discussed  in  section  2.9.  These  functions  play  an  im¬ 
portant  role  in  photomechanics  and  will  be  given  the 
names  : 

D(t)  -  optical  creep  function  , 

E ( t )  -  inverse  optical  creep  function. 

4.2.  Stress  and  Strain  Birefringence. 

Consider  the  possibility 

B(t)  =  Bq  h(t)  ,  (9) 
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where  h(t)  is  the  unit  step  function.  Then 


R.  .  =  2B 

1  J  o 

e .  . 

iJ  > 

(10) 

E(t)  =  | 

o 

G(t)  , 

(11) 

D(t)  =  B 
v  J  o 

J(t)  . 

(12) 

In  this  case  the  deviatoric  part  of  the  dielectric 
tensor,  and  therefore  the  observed  fringe  order  in  a 
polariscope,  is  directly  related  to  the  strain.  A 
material  for  which  (9)  holds  may  be  termed  s tr ain- 
jbi  ref ringent .  The  isoclinic  angle  then  coincides 
with  the  principal  axis  of  strain. 

Consider  the  possibility 

D(t)  =  Do  h(t)  .  (13) 

Then, 


R.  .  =  D  S.  .  , 

1J  o  1J  ' 

(14) 

B(t)  =  Dq  G(t)  , 

(15) 

E(t)  =  h(t)  . 

(16) 

o 


Such  a  material  may  be  called  stress-birefringent. 

The  fringe  order  then  gives  the  difference  in  prin¬ 
cipal  stresses  and  the  isoclinic  angle  coincides  with 
the  principal  axes  of  stress. 


Consider  an  elastic  material  which  is  stress- 
hiref ringent .  Combining  the  relations  of  section  2.2 
with  (13)  -  (16)  leads  to  a  relation  of  the  form  (10). 
Thus,  a  stress-biref ringent  elastic  material  is  also 
strain-birefringent.  Similarly,  a  strain-biref ringent 
elastic  material  is  also  stress -biref ringent . 

Consider  the  possibility 

D(t)  =  a  h(t)  +  1  b  J(t)„  (17) 


where  a  and  b  are  constants.  Then, 


R..  =  a  s..  +  b  e.. 
ij  ij  iJ 


(18) 


Thus,  the  general  expression  includes  the  possibility 
that  the  fringe  order  depends  on  a  combination  of  the 
stress  and  strain.  The  relation  (18)  is  often  in 


fair  agreement  with  experiment. 

The  list  of  such  special  cases  is  endless.  The 
relations  obtained  may  describe  some  materials  accurate¬ 
ly  enough  for  some  applications,  but  they  are  certainly 
not  universal.  The  relation  (18)  is  particularly  useful 
in  qualitative  reasoning,  i.e.,  to  exercises  of 


"physical  intuition". 


4.3.  Slow  on  Rapid  Motion. 

The  constitutive  relations  (1)  are  of  the  same 


form  as  the  stress  strain  laws  of  viscoelastic 


55 


materials.  Thus,  the  same  approximation  theorems  hold 
for  both  relations. 

For  sufficiently  rapid  straining,  by  (3), 

E..  =  2  B(0)  e..  .  (19) 

The  material  is  strain  biref ringent .  Since,  by  (57)  of 
chapter  2,  the  material  is  elastic  under  such  straining, 
the  material  is  also  stress-biref ringent . 

For  slow  straining,  by  (3), 

R.  .  =  2  B(eo)  e.  .  .  (20) 

ij  ij 

The  material  is  strain-biref ringent.  Since,  by  (59) 
of  chapter  2,  the  material  is  elastic  for  slow  straining, 
it  is  also  stress-biref ringent . 


4.4.  Temperature  Dependence. 

Each  of  the  material  properties  functions  of 
section  4.1  is  temperature  dependent.  From  the  ex¬ 
perimental  observations  of  the  similarity  of  mechanical 
and  optical  creep,  it  may  be  expected  that  materials 
which  are  thermorheologically  simple  in  their  mechani¬ 
cal  behavior  are  described  by  similar  relations  for  the 
optical  behavior. 

Let  the  reduced  time  be  defined  by 
C (t)=  b(T(t) )  dr  . 


The  function  b(T)  of  temperature  may  be  termed  the 
optical  shift  factor.  An  optically  and  thermo-rheologi- 
cally  sijnple  material  is  then  defined  by  the  constitu- 
tive  relation: 

R.  .  =  2B(C(t)-C(T))e  .(t)  dr.  (22) 

1  J  O  -L  J 

This  generalizes  the  relation(3). 

It  may  happen  that  the  mechanical  and  optical 
shift  factors  are  the  same.  However,  this  must  be 
determined  experimentally. 
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5.  BASIC  FORMULAE  FOR  PLANE  STRESS. 

5.1.  Normal  Incidence. 

Consider  a  slice  of  material  of  thickness  h  which 
experiences  a  state  of  plane  stress  described  in 
section  2.4.  The  state  of  stress  is  to  be  determined 
by  observing  the  isoclinic  and  isochromatic  fringe 
patterns  in  a  polariscope. 

The  fringe  order  n  and  the  isoclinic  angle  6  are 

observed  as  a  function  of  time.  The  results  of  chapter 

3  are  used  to  derive  information  about  the  dielectric 

tensor.  Although  the  solution  obtained  there  is  exact 

only  when  K.  .  is  constant,  it  will  be  used  here  as  an 

approximate  solution  when  K.  .  varies  in  both  time  and 

r  J 

space.  This  approximation  will  be  satisfactory  as 
long  as  the  time  scale  of  the  mechanical  behavior  is 
long  compared  to  the  time  ,  about  10  ^  seconds; 
that  is,  the  time  required  for  light  to  traverse  the 
model.  The  results  of  chapter  4  are  then  used  to 
derive  information  about  the  stress  tensor. 

The  x-y  axis  may  be  chosen  arbitrarily.  The 
principal  axes  of  stress  are  then  located  at  an  angle  5 


from  the  x  toward  the  y  axis: 
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tan  2?  = 


xy 
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CTx  "  CTy 


xy 
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Let  Q  denote  the  angular  orientation  of  the 


principal  axis  of  K.  .  in  the  plane.  Then, 

u 


R  -  R  =  K 
x  y  x 


Ky  =  (Kx  -  K2)  cos  2Q 


R  =  K  =  |(K,-  K0)  sin  2Q 
xy  xy  1  2 


By  (42)  of  section  3.5, 
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N  c 

V  Ry  =  <pr  P2>  003  2«  - 
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R— -  =  I  TT—  (Pi-Pp5  Sin  2C  ‘ 


xy 
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By  (51)  and  (52)  of  section  3.5, 


4ttN  c 

R  -  R  =  - r -  n  cos  2Q  , 

x  y  uuh  ^ 


R 


2rrN  c 
o 


xy  uuh 


n  sin  2Q  , 


and  n  is  the  iso chromati c  fringe  order  . 


By  (5)  of  section  4.1, 

n  cos  2Q  =  Jq  Y  ( t-  t  )  C(o1(t)-  ct2(t))cos 

n  sin  2Q  =  Y ( t- t )  [ct-^t)-  ct2(t))  sin 


where  , 


Y  = 


tuh  3 


4rrNoc 


,  (11) 
(12) 

(13) 

(14) 

(15) 

(16) 

2§(T)]dT,(17) 
2? (t) ]dT; (18) 

(19) 
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By  (8)  of  section  4.1> 


(a]_ -a2)  cos  2?  =  J^(t-T)  ^  [n(T)cos  2C(t)]  dT,  (20) 
(a1-a2)  sin  2?  =  §(t-T)  [n(T)  sin  2C(t)]  dT;(21) 

wh  ere , 

4ttN  c 

§  =  - 2 -  E  „  (22) 

u)h 

Equations  (20)  and  (21)  can  be  used  to  determine 
the  difference  in  principal  stresses  and  the  orientation 
of  principal  axes  of  stress  at  a  point  from  the  observed 
history  of  fringe  order  and  isoclinic  angle  at  the 
point  with  normal  incidence  of  light. 

The  principal  axes  of  stress  and  strain  coincide 
with  the  isoclinic  angle  if  and  only  if  the  isoclinic 
angle  is  constant.  Otherwise  all  three  sets  of  prin¬ 
cipal  axes  are  different. 

When  the  isoclinic  angle  is  independent  of  time, 
the  equations  (20)  and  (21)  give  the  result: 

al"a2  =  Jo  n(T)  dT  *  ( 2 3) 

Invers  ely 

n  =  Jo  dT  [CTl(T)‘  a2 ( T ) ~^dT  •  (24) 
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5.2.  Oblique  Incidence. 

The  observation  of  the  history  of  the  fringe 
pattern  with  the  incident  light  directed  normal  to  the 
slice  gives  the  difference  in  principal  stresses  and 
the  orientation  of  the  axis  of  principal  stress.  The 
values  of  the  two  principal  stresses  are,  therefore, 
not  determined  solely  by  the  observation  of  fringe 
order  and  isoclinic  angle  at  a  point  with  normal  inci¬ 
dence. 

If  the  principal  axes  of  stress  are  found  to  have 
the  same  orientation  at  all  times,  then  the  method  of 
oblique  incidence  can  be  conveniently  used  to  give  a 
third  independent  piece  of  information:  From  the  com¬ 
bined  observations  with  normal  and  oblique  incidence 
the  two  principal  stresses  can  be  calculated.  If  the 
principal  axes  of  stress  are  not  independent  of  time 
the  method  of  oblique  incidence  can  still  be  used  but 
the  relations  are  more  complicated  and  will  not  be  given 
in  this  paper. 

The  z-axis  is  normal  to  the  slice  and  the  x,y- 

axes  are  taken  along  the  principal  axes  of  stress.  The 

case  when  the  principal  axes  of  stress  do  not  change 

with  time  will  be  considered.  By  (5)  of  section  4.1, 

the  principal  axis  of  K.  .  do  not  rotate.  Thus,  the  x, 

i  j 

y-axes  are  principal  axes  of  the  dielectric  at  all 
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times.  As  an  example,  suppose  that  the  incident  light 
is  directed  normal  to  the  y-axis  and  making  an  angle  a 
with  the  z-axis.  The  fringe  order  can  be  related  to 
the  dielectric  tensor  by  (72)  -  (74)  of  section  (3). 

The  dielectric  tensor  is  related  to  the  stress  tensor 
by  (4)  and  (5)  of  section  4.  Combining  these  relations 
for  the  case  when  the  x,y-axes  are  the  principal  axes 
of  both  stress  and  dielectric  tensors  leads  to  the 
relations  0 

n=  Jo  'Kt“T)  [^(t)  cos2e  -  o  dT  •  (25) 

cos  3 
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6.  PREPARATION  OF  MODEL  MATERIAL. 

It  is  usually  desired  to  determine  the  stress  in  a 
viscoelastic  solid  from  observations  on  a  geometrically 
similar  model.  The  experimental  procedure  determines  the 
stresses  in  the  model.  The  correlation  with  the  prototype 
generally  requires  that  the  viscoelastic  properties  of  the 
model  be  similar  to  those  of  the  prototype:  The  Poisson's 
ratio  in  creep  v(t)  must  be  the  same  and  the  tensile  creep 
compliances  must  be  proportional  for  the  model  and  prototype 
materials.  Differences  in  Poisson's  ratio  do  not  generally 
lead  to  significant  errors.  The  compliances  are  always 
proportional  for  elastic  materials;  and  so,  the  choice  of 
model  material  does  not  present  a  problem  in  photoelasticity. 

It  is  generally  necessary  to  manufacture  a  new  material 
for  each  viscoelastic  experiment.  One  such  material  is 
described  in  this  section.  Many  other  materials  could  be 
used. 

An  epoxy  mixture  made  from  components  manufactured  by 
the  CIBA  Co.  can  be  made  to  give  considerable  variation 
in  viscoelastic  properties.  The  material  described  here  is 
made  of  CIBA  Araldite  502,508,  and  963,  in  the  proportions 
45/45/10.  Two  different  batches  referred  to  as  A  and  B, 
having  the  same  formula,  are  described. 

A  mold  is  prepared  to  cast  sheets  of  the  material  /  inch 
thick  and  two  square  feet  in  area:  Sheets  of  plate  glass  are 


64 


used  for  the  two  faces  of  the  mold,  £  inch  surgical  tubing 
forms  a  gasket  between  the  sheets  on  three  sides,  the  thick¬ 
ness  is  controlled  by  x  inch  spacers  placed  between  the  glass 
faces  and  outside  of  the  gasket. 

This  epoxy  formulation  is  a  very  tenacious  adhesive 
and  care  must  be  taken  to  insure  mold  removal:  The  surfaces 
must  first  be  cleaned  with  a  solvent.  Three  light  coats 
of  a  Carnuba-base  wax  are  then  applied,  allowing  each  coat 
to  partially  dry  and  polishing  it  before  application  of  the 
next  coat.  Finally  all  surfaces  are  sprayed  with  a  polyvinyl 
alcohol  solution  which  when  properly  applied  and  dried  forms 
a  smooth  water  soluble  film.  Some  experience  is  required  to 
evolve  a  spraying  technique  which  gives  a  glassy  smooth 
surface.  The  Carnuba  base  wax  and  the  polyvinyl  alcohol 
solutions  are  available  through  most  commercial  fiberglass 
supply  houses . 

At  room  temperature  the  viscosity  of  the  uncatalyzed 
resin  components  is  such  that  small  air  bubbles  will  become 
trapped  in  the  mixture  while  stirring,  resulting  in  a  poor 
casting.  If  the  mixture  is  placed  in  a  vacuum  some  bubbles 
will  be  drawn  out  but  the  mixture  will  catalyze  before  all  of 
the  bubbles  are  gone.  However,  if  the  resin  components  are 
heated  to  130°F  prior  to  mixing  the  viscosity  is  lowered 
sufficiently  to  allow  entrained  air  to  float  off  before 
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solidification.  Keating  the  components  reduces  the  pot  life 
of  the  mixture  to  less  than  five  minutes  after  the  hardener 
is  added.  Final  mixing  and  pouring  into  the  mold  must  be 
done  quickly.  A  paddle  mounted  at  the  end  of  a  rod  and 
driven  by  a  slow  speed  drill  motor  facilitates  mixing. 

The  resin  is  cured  in  an  oven  for  12  hours  while  the 
temperature  is  maintained  at  160  ±  5°F.  The  mold  is  taken 
from  the  oven  and  the  clamps  removed.  The  rubber- tube 
gasket  is  stripped  off  and  a  razor  blade  is  used  to  cut 
around  the  edges  between  the  cast  sheet  and  the  glass.  Wooden 
wedges  are  slowly  pressed  between  the  edges  of  the  glass 
faces  while  the  sheet  is  held  under  running  water. 

Careful  prying  will  remove  the  mold  after  about  five 
or  ten  minutes.  The  newly  cast  sheet  is  washed  to  remove  the 
mold  release  agent  and  placed  on  a  flat  paper  or  teflon 
coated  surface  to  relax  to  equilibrium. 

After  the  sheet  has  relaxed  for  a  day,  models  and  speci¬ 
mens  may  be  laid  out  on  the  sheet.  These  shapes  are  cut  from 
the  sheet  on  a  band  saw  at  high  speed  with  a  skip- tooth  blade. 

Machining  accurate  models  from  this  material  presents 
some  problems.  If  the  material  is  clamped  lightly  in  a  vice 
or  a  chuck  it  relaxes  and  becomes  loose.  If  the  vice  or 
chuck  is  tightened  enough  to  hold  the  part,  the  part  becomes 
distorted  so  that  a  flat  surface  or  a  square  edge  after 
removal  from  the  jig  is  no  longer  flat  or  square.  Another 
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difficulty  appears  when  a  heavy  cut  is  made:  The  tool  mashes 
the  material  under  it  slightly  so  that  after  relaxation  the 
part  is  larger  than  expected. 

These  difficulties  have  "been  overcome  hy  using  double- 
sided  adhesive  tape  to  hold  the  model  to  the  face  plate  for 
turning  operations  and  by  taping  the  jaws  and  bottom  of  the 
vice  for  milling  operations.  If  sharp  corners  are  required, 
the  model  may  be  taped  on  both  faces  and  sandwiched  between 
1/8  inch  plexiglas,  the  sandwich  being  then  taped  on  both 
edges  and  very  lightly  clamped  in  a  vice.  If  a  particular 
dimension  is  desired,  it  should  be  approached  slowly  with 
fine  cuts.  When  the  machine  has  reached  the  final  setting, 
the  specimen  is  allowed  to  relax  for  15  to  30  minutes  and  the 
last  cut  is  made  at  the  same  machine  setting.  This  process 
can  be  repeated  if  extreme  accuracy  is  required. 

Flat  surfaces  and  edges  are  cut  with  a  sharp  fly  cutter 
running  at  high  speed.  Circular  section  can  be  turned  in  a 
lathe  at  high  speed  using  a  sharp  pointed  tool  with  about 
35°  rake.  Holes  may  be  drilled  if  the  model  is  sandwiched 
between  plexiglas  and  the  drill  is  sharp  and  turned  at  high 
speed . 

For  some  models  it  is  helpful  to  cool  the  model  in  a 
refrigerator  immediately  before  machining.  The  material  is 
quite  stiff  at  40°F. 

The  models  have  been  machined  dry.  If  a  coolant  is  to  be 
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used,  experiments  should  he  made  to  determine  if  the  coolant 
accelerates  the  growth  of  edge  fringes.  Water,  for  example, 
will  considerably  accelerate  the  growth  of  edge  fringes. 

Finished  models  stored  at  normal  room  conditions  will 
be  useable  for  several  weeks  or  more  before  edge  fringes  be¬ 
come  excessive.  Carnuba  wax,  acetone,  n-heptane,  methyl- 
enedichloride  and  water  will  accelerate  edge  fringe  growth. 

A  dessicating  atmosphere  will  retard  edge  fringe  growth. 

The  photoviscoelastic  material  is  calibrated  in  a  ten¬ 
sile  creep  test.  A  typical  tensile  specimen  thick  by 
0.6"  long  is  shown  in  Fig.  1.  Usually  six  or  eight  cali¬ 
bration  specimens  are  cut  from  each  sheet  of  model  material. 
Experiments  have  shown  that  specimens  cut  from  different 
locations  and  orientations  in  the  sheet  yield  the  same 
creep  behavior. 

The  entire  batch  of  specimens  may  be  stacked  on  edge 
and  milled  to  a  uniform  width  in  one  set-up  with  a  fly  cutter. 
Aluminum  tabs  are  bonded  to  each  end  of  each  specimen  with 
Eastman  910  and  the  specimens  are  placed  on  edge  on  a  flat 
surface  to  relax  before  testing.  By  having  six  or  eight 
specimens  any  one  specimen  has  enough  time  to  relax  between 
tests  while  the  other  specimens  are  tested. 
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7 


MATERIAL  CHARACTERI ZATION 


7.1.  Creep  Test. 

The  mechanical  and  optical  properties  are  deter¬ 
mined  by  the  tensile  creep  test;  A  constant  uniaxial 
stress  is  applied  and  the  strain  and  birefringence 
are  recorded,  figure  2.  The  load  is  applied  to  the 
specimen  by  releasing  a  weight  and  total  change  in 
length  is  measured  with  a  differential  transformer. 

The  load  is  measured  with  a  beam  instrumented  with 
strain  gages.  The  time  interval  during  loading  is 
approximately  0.2  seconds.  A  calibrated  clip  gage, 
made  by  bending  a  0.010  inch  thick  strip  of  aluminum 
and  instrumenting  it  with  strain  gages,  is  used  to 
measure  the  changing  width  of  the  specimen  as  it  creeps. 

All  strain  gages  and  differential  transformers 
are  powered  by,  and  their  outputs  recorded  on,  a 
6-channel  Brush  recorder  with  a  calibrated  time  base. 

The  fringe  growth  is  also  automatically  recorded: 
The  output  of  a  mercury-arc  lamp  is  filtered  with  a 
Wratter  77A  filter  which  passes  the  green  line 
(5460  A)  of  mercury.  This  light  is  passed  through  a 
polarizer  (Polaroid)^ through  the  creep  test  specimen, 
through  a  second  polarizer  with  axis  at  right  angles 
to  the  first  and  then  into  a  photomultiplier  tube. 

The  output  of  the  photomultiplier  tube  shows  the  periodic 
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increasing  and  decreasing  of  light  intensity.  The 
maximums  and  minimums  of  this  curve  are  well  defined 
and  represent  increasing  integral  fringe  orders 0 

The  model  is  enclosed  in  a  box  which  has  electric 
resistance  heaters,  automatically  controlled  to  maintain 
a  set  temperature  ±  0.2°F  over  the  range  from  room 
temperature  to  150°F.  Circulating  fans  within  the  test 
chambers  insure  a  uniform  temperature. 

The  glass  transition  temperature  for  the  model 
material  is  in  the  region  of  room  temperature.  Small 
changes  in  temperature  (5°F)  cause  marked  changes  in 
mechanical  behavior.  At  130  to  150°F  the  relaxation 
time  is  very  short  and  the  material  is  rubbery  while 
at  32°F  it  is  quite  rigid. 

For  the  tensile  creep  test,  the  principal  axes  of 
stress,  strain,  and  refraction  coincide.  Let  a  denote 
the  tensile  stress  .  Equation  (24)  of  section  5.1  leads 
to 

+  (t)  =  .  (i) 


From  (12)  of  section  2.1, 
e  (t) 

n(t)  -  -j-  . 

From  (13)  of  section  2.1, 
e  (t) 

>»(*)  *  -ritT  • 

X  ' 


(2) 

(3) 
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Then  J(t)  and  B(t)  can  he  calculated  from  (16)  and  (17) 
of  section  2.1. 

J  =  2(1+  v)  D,  (4) 

B  =  3(1-  2v)  D.  (5) 

The  constant  stress  a  is  taken  to  he  the  load 
divided  hy  the  undeformed  cross-section  area.  The 
successive  peaks  on  the  light  intensity  curve  corres¬ 
pond  to  integer  values  of  the  fringe  order  n.  By 
observing  the  time  of  occurrence,  the  optical  creep 
function  t|r(t)  can  be  determined  from  (l).  A  graph  of 
\|[  determined  at  room  temperature  for  mix  A  is  shown  in 
Fig.  3.  The  results  of  tests  at  various  temperatures 
are  shown  in  Fig.  4.  for  mix  B.  These  two  mixes  were 
of  the  same  proportions  of  raw  material  hut  show  a 
slight  difference  in  their  mechanical  and  optical 
properties. 

The  longitudinal  strain  was  taken  to  he  the  total 
change  in  length  divided  by  the  original  length.  The 
tensile  creep  compliance  D(t)  computed  by  (2)  is  shown 
in  Fig.  5,  and  the  transverse  contraction  ratio  v(t)  is 
shown  in  Fig.  6.  Because  of  the  small  change  in 
thickness  the  lateral  strains  are  difficult  to  measure, 
and  the  large  band  of  experimental  data  in  Fig.  6  is 
due  to  this  experimental  error  rather  than  any  real 
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variation  of  material  properties.  The  shear  creep 
compliance  J(t)  and  hulk  creep  compliance  B(t)  cal¬ 
culated  from  (4)  and  (5)  are  shown  in  Fig.  7  and  8. 

In  the  tests  at  higher  temperatures  of  mix  B,  only  the 
longitudinal  strain  was  measured.  The  tensile  creep 
compliance  is  shown  in  Fig.  9. 

If  the  dependence  of  the  final  data  on  temperature 
is  neglected  in  Fig.  4,  then  the  temperature  dependence 
of  the  creep  functions  over  this  small  change  of  tem¬ 
perature  can  he  reasonably  well  described  by  a  tem¬ 
perature  shift  factor.  The  optical  and  mechanical  shift 
factors  are  shown  in  Fig.  10. 

The  optical  creep  function  is  also  reasonably 
well  described  by  the  empirical  relation  (18)  of  section 
4.2  with  a  =  18.4  x  10  ^  in/lb^  and  b  =36.8. 

The  resulting  error  is  20%  for  some  ordinates. 

In  order  to  analyze  the  photoviscoelastic  data  it 
is  necessary  to  calculate  the  inverse  of  the  optical 
creep  function  \|f(t)  related  to  $(t)  by 

Jo  $(t)  t(t-T)  dT  =  t  . 

This  equation  is  solved  by  a  numerical  technique  in 
conjunction  with  the  digital  computer  as  described  in 
section  2.9.  The  results  are  shown  in  Fig.  11  for  one 
cas  e . 
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7.2.  Thermal  Properties . 

The  coefficient  of  thermal  expansion  was  measured 
in  the  jig  shown  in  Figure  12.  The  specimen  rests  on 
a  Teflon- coat ed  surface  of  the  jig  and  one  end  is  bonded 
to  the  jig.  A  metal  tab  is  bonded  to  the  free  end  of 
the  specimen  and  a  fine  wire  soldered  to  the  tab  leads 
to  one  terminal  of  an  ohmmeter.  The  ohmmeter  registers 
a  short  circuit  when  the  micrometer  is  brought  to  touch 
the  free  end  of  the  specimen.  By  using  this  system  the 
micrometer  measurements  may  be  repeated  to  within 
±  0.0002.  Such  accuracy  cannot  be  obtained  by  mechani¬ 
cal  feel  due  to  the  low  modulus  of  the  model  material. 

The  jig  was  placed  in  an  oven,  allowed  to  come  to 

equilibrium  at  several  temperatures  and  changes  of 

length  were  measured.  The  hysteresis  was  determined 

by  approaching  a  given  temperature  from  above  and  below 

and  was  small.  The  temperature  was  varied  between 

70  and  110°F  and  the  coefficient  of  thermal  expansion 

was  calculated  to  be  1.06  x  10  ^  - — -  ±  6%  . 

o  _ 

m  -  F 

Correction  was  made  for  the  expansion  of  the  measuring 

jig. 

The  coefficient  of  thermal  diffusivity  was  cal¬ 
culated  by  measuring  the  temperature  history  at  the 
centers  of  two  slugs  immersed  in  a  constant  temperature 
bath.  The  temperatures  were  measured  with  iron- 
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constantan  thermocouples  cast  into  the  slugs.  The 
constant  temperature  bath  was  water  in  a  large 
laboratory  thermos  bottle. 

The  slugs  and  the  bath  were  allowed  to  reach  their 
separate  equilibrium  temperatures  as  determined  by  the 
thermocouple  readings.  The  slugs  were  then  immersed 
in  the  bath  and  stirred  vigorously.  The  temperature 
history  being  recorded  on  a  calibrated,  self-balancing 
potentiometer.  A  rectangular  and  a  cylindrical  slug 
were  each  tested  several  times  to  insure  reproduci¬ 
bility.  The  thermal  diffusivity  was  calculated  for 
both  geometries  at  several  points  in  the  temperature 
history  and  was  32.9  x  10  ^  ft'Vhour  ± 

The  specific  heat  was  determined  by  a  calorimeter 
test  and  was  0.503  BTU/lb°F  ±  8%.  The  density  of  the 
model  material  is  O.O46  lb/cubic  inch.  From  these 
values  the  thermal  conductivity  was  computed  to  he 
0.131  BTU/hr  ft°F  ±  10%. 
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SHEAR  TEST 


B.lo  Specimen  and  Apparatus, 

A  sheet  was  bonded  to  metal  blocks  along  opposite 
edges  as  shown  in  Fig.  13.  A  state  of  nearly  plane 
stress  is  obtained  in  the  sheet  by  loading  the  blocks. 

If  the  load  on  the  block  is  constant,  the  problem  is  one 
of  proportional  loading  as  described  in  section  2.3. 

The  stresses  will  be  constant  if  the  time  variation  of 
Poisson's  ratio  can  be  neglected.  Experience  with 
elasticity  problems  indicates  that  the  observed  varia¬ 
tion  of  v  should  not  have  a  significant  effect  on  the 
stress  state.  Therefore  this  configuration  was  chosen 
for  some  preliminary  tests. 

The  polariscope  used  is  shown  in  Fig.  14.  It  is 
a  versatile  piece  of  equipment  for  general  use  in  photo¬ 
viscoelasticity.  It  differs  from  the  ordinary  photo¬ 
elastic  bench  in  that  the  polari zer- analyzer  combination 
can  be  rapidly  rotated  by  electric  motors. 

Three  ^Light  sources  may  be  used:  a  point  source 
(d.c.  mercury  arc  with  a  Wratten  77A  filter),  a  distri¬ 
buted  source  (a.c.  mercury  arc  and  a  green  fluores  cent),  and 
diffused  source.  Condensing  lenses  are  used  with  the 
point  source  light  to  get  parallel  beam  through  the 
specimen.  The  lenses  are  removed  from  the  light  path 
when  the  distributed  light  sources  are  used. 
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Each  lens  stand  is  fitted  with  a  ball-hearing 
carrier  to  support  the  rotating  polaroid  around  its 
outside  edge.  The  polaroid  is  driven  by  a  stepping 
solenoid  through  a  positive  drive,  link  chain.  When 
the  solenoid  is  pulsed  the  polaroid  advances  9°.  The 
pulse  rate  is  continuously  adjustable  from  5  pulses 
per  second  to  1  pulse  per  2  minutes. 

A  16  mm  Cine  Special  camera  is  used  to  record  the 
fringe  patterns.  Mounted  next  to  the  model  and  in  the 
field  of  view  of  the  camera  is  a  digital  timer  with 
l/lO  second  minimum  reading.  Incorporated  in  the  timer 
is  a  numbered  wheel  driven  by  a  small  rotating  solenoid 
which  is  tied  to  the  driver  solenoids  and  indicating  the 
angle  of  the  polarizing  axes. 

During  the  first  minute  after  loading  the  camera 
is  run  at  16  frames  per  second  to  capture  the  rapidly 
changing  fringe  pattern.  The  frame  rate  is  then  re¬ 
duced  to  1  frame  per  second.  As  the  fringe  patterns 
approach  equilibrium  the  frame  rate  is  correspondingly 
reduc  ed . 

A  model  is  mounted  in  grips  in  the  loading  frame 
and  the  upper  grips  are  connected  to  a  load  cell  through 
a  universal  joint.  Two  load  cells  of  differing  sensi¬ 
tivities  are  available.  Four  loading  patterns  may  be 
used:  constant  load,  constant  strain,  constant  strain 
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rate  and  sinusoidal  strain. 

Constant  loads  are  applied  by  hanging  weights 
on  the  specimens  and  constant  strains  are  applied  by 
using  a  heavy  weight  to  force  the  loading  mechanism 
against  a  present  stop.  Constant  strain  rate  is  applied 
to  the  model  by  driving  the  lower  grips  with  a  220  volt, 
3  phase  motor  operating  through  a  gear  change  set.  By 
changing  the  gears  the  strain  rate  may  be  varied  between 
0.0413  and  5.75  inches/minute.  A  spring  loaded  can  and 
tappet  mechanism  supplies  the  sinusoidal  strains. 

Figure  A. 6. 2.  The  mechanism  is  driven  by  an  electric 
motor  through  a  Zeromax  speed  reducer.  The  frequency 
is  continuously  variable  from  0  to  400  cycles/minute 
and  the  maximum  amplitude  may  be  varied  between  0  and 
0.2  inches,  by  inserting  appropriate  cans.  All  of  the 
loading  mechanisms  are  mounted  in  a  rigid  frame. 

Strains  are  measured  with  differential  transformers 
and  read  out  on  a  Brush  recorder  or  an  oscilloscope 
equipped  with  a  camera. 

8.2.  Creep  in  Shear. 

In  the  first  test,  the  model  was  "sheared"  by 
constant  loads  As  shown  in  Section  2.3,  if  the 

variation  in  v  is  neglected,  the  stress  will  be  constant 
and  be  the  same  as  if  the  material  were  elastic.  The 
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stress  state  in  an  elastic  material  was  determined  by 
standard  photoelastic  means.  It  was  found  that  the 
shear  stress  varied  along  the  centerline  as  shown  in 
Fig.  15.  By  symmetry,  the  point  in  the  center  of  the 
specimen  experiences  a  pure  shear  stress,  i.e.  cx=6y=  0. 
Thus,  constant  loads  P1  are  expected  to  give  a  constant 
state  of  stress 
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The  particle  at  the  center  experiences 
By  the  formula  of  Section  5.1,  the 
of  fringe  order  is 


n ( t )  =  2t  (t)  . 
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are  shown  in  Fig.  16.  They 
value.  The  small  deviation 
the  large  shear  strains  which 
linear  behavior. 


8.3.  Rotation  of  Principal  Axes. 

It  was  then  decided  to  demonstrate  the  fact  that 
the  isoclinic  angle  does  not  coincide  with  the  principal 
axes  of  stress.  By  first  applying  a  load  ?2  and  then. 


78 


the  prin- 


at  a  later  time  t  ,  applying  the  load  P^, 
cipal  axes  of  stress  can  he  made  to  rotate.  The  iso¬ 
clinic  angle  can  he  simultaneously  measured. 

The  weight  of  the  frame  alone  caused  some  shear 
stress  which  may  he  regarded  as  having  reached  an 
equilibrium  state.  By  (5)  of  section  4.1, 
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where  tq  =  0.30  psi. 

The  loads  result  in  a  state  of  stress  at  the 

center  with  t  =0.  The  stresses  a  and  cr  determined 
xy  x  y 

hy  photoelasticity  using  oblique  incidence  were: 


Thus 


°x  -  X-95  p2' 

a  ~  0.22  Pp , 
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t  =  o . 
xy 
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The  shear  load  applied  at  time  tQ  results  in 
stresses  : 


xy 


=  0  , 
y  -  0  • 


(6) 


From  the  photoelastic  test  t  =  2.5  P^. 
Thus  , 
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xy 


D(t-to)  T  , 


(7) 


Since  the  problem  is  linear,  the  response  to  the 
three  stress  states  can  be  superimposed.  The  combined 
values  of  the  tensor  R^.  for  t  >  tQ  are  determined 
by  adding  (3),  (5),  and  (7): 


R  =  D(eo)  t  +  D(t-t  )t 
nxy  O  O 

R  -  R  D  ( t )  o 
x  y 

The  isoclinic  angle  is  given  by 

2R 

tan  2 C  =  p  _r 

x  y 


(8) 
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Thus  , 
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(10) 


The  angle  C  predicted  by  (10)  is  compared  with  the 
observed  isoclinic  angle  in  Fig.  17.  The  agreement 
is  within  experimental  error. 
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CIRCULAR  DISK 


9.1.  Specimen  and  Apparatus. 

Typical  models  are  shown  in  Fig. 18  and  Fig.  19. 
The  outer  boundary  is  circular  but  the  inner  boundary 
may  be  an  irregular  shape.  The  model  may  be  heated 
along  the  inner  boundary  and  loaded  by  a  uniform 
pressure  along  the  outer  boundary. 

The  loading  jig  is  pictured  in  Fig.  20.  Pressure 
is  applied  through  a  0.010"  thick  latex  diaphragm. 

The  diaphragm  is  designed  to  have  enough  slack  to 
follow  the  model  as  it  deforms  without  stretching  thus 
maintaining  uniform  pressure. 

The  uniformity  of  the  loading  and  the  pressure 
efficiency  of  the  diaphragm  were  determined  with  the 
elastic,  Hysol  4485  specimens  pictured  in  Fig.  19. 
Knowing  the  elasticity  solution  for  the  circular  disk 
and  the  photoelastic  fringe  constant  for  the  calibra¬ 
tion  model  material  one  may  calculate  the  effective 
pressure  acting  on  the  model.  This  compared  to  the 
gage  reading  of  the  air  pressure  acting  on  the  dia¬ 
phragm  establishes  the  pressure  efficiency  of  the  jig. 
Non- uniformity  in  the  pressure  around  the  edges  of  the 
disk  will  cause  the  fringe  pattern  to  depart  from  con¬ 
centricity.  A  measure  of  this  non- uniformity  is  the 
per  cent  variation  in  fringe  order  around  an  imaginary 
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concentric  circle  drawn  on  the  model.  A  picture  of 
a  loaded  calibration  specimen  is  shown  in  Fig.  21. 

The  pressure  is  uniform  to  within  10%  and  the  jig 
efficiency  (effective  pressure/gage  pressure)  is  68%. 

The  model  is  sandwiched  in  the  jig  between  two 
Plexiglas  guards.  Polaroid  is  bonded  to  the  model 
side  of  each  guard.  These  guards  are  not  loaded. 

They  serve  to  contain  the  thin  latex  diaphragm,  to 
prevent  the  low  modulus  model  from  buckling,  and  to 
insulate  the  heated  model.  A  copious  layer  of  thin 
oil  separates  the  model  from  the  guards.  Calibration 
tests  show  that  the  friction  is  insignificant. 

The  model  is  heated  along  the  inner  edge  by  a  re¬ 
sistance  wire  heater  fitted  into  the  opening.  In  the 
tests  conducted,  the  model  was  heated  to  an  equilibrium 
temperature  and  a  preliminary  test  was  run  to  determing 
the  temperature  dis  tribution.  A  dummy  model  and  guard  of 
the  same  materials  and  dimensions  as  the  test  model  and 

guard  were  placed  in  the  loading  jig.  Holes  0.020"  in 

\ 

diameter  were  drilled  through  the  guard  and  into  the 
model  at  several  stations.  These  holes  were  filled  with 
oil  to  minimize  their  effect  on  the  temperature  dis¬ 
tribution.  The  model  was  heated  to  equilibrium  with 
constant  power  input  to  the  heater.  A  small  iron- 
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constanan  thermocouple  probe  was  placed  sequentially 
into  the  holes  to  measure  the  temperature.  Several 
tests  showed  that  the  temperature  distribution  could 
be  repeated  to  within  ±  1°F.  The  model  and  apparatus 
used  for  this  test  is  shown  in  Fig.  22. 

9 .. 2 .  Exterior  Pressure . 

The  model  of  Fig.  18  was  loaded  with  a  step 
input  (approximately  0.1  second  rise  time)  of  pressure 
which  was  maintained  constant.  The  changing  fringe 
patterns  were  recorded  on  a  16  mm  movie  camera  running 
at  16  frames  per  second,  shown  in  Fig.  23. 

It  was  shown  in  Section  2.4  that  the  stress  field 
in  a  plane  stress  problem  with  stress  boundary  conditions 
and  uniform  temperature  is  the  same  as  if  the  material 
were  elastic.  In  particular,  the  stress  is  independent 
of  the  material  properties.  Thus,  if  the  unheated 
model  is  subjected  to  constant  external  pressure  8,  the 
stress  at  any  point  should  be  constant:  Each  particle 
experiences  creep  at  constant  stress.  The  isoclinic 
angle  at  each  point  is  then  constant  and  coincides  with 
the  principal  axes  of  stress.  With  the  coordinate  axes 
along  the  principal  axes,  equation  (24)  of  section  5.1 
leads  to 

n(t)  =  ilf(t)  [a1  -  a2l  .  (!) 


84 


Since  \|i(t)  is  known  from  the  tensile  test,  the  expected 
fringe  order  can  he  calculated,  to  within  a  constant 
factor. 

The  actual  fringe  pattern  was  recorded  with  a 
movie  camera.  A  typical  frame  is  shown  in  Fig.  24. 

The  expected  fringe  order  is  compared  with  the  observed 
fringe  order  at  several  points  in  Fig.  25,  and  the 
agreement  is  satisfactory.  The  stress  is  constant 
and  the  stress  at  the  star  root,  point  A,  is  found 
to  be  41  psi  or  3.9  p.  This  agrees  with  the  photo¬ 
elasticity  solution. 

9.3.  Exterior  Pressure  with  Temperature  Gradient. 

The  model  was  heated  until  a  steady  non-uniform 
temperature  distribution  was  reached.  This  temperature 
was  maintained  throughout  the  loading.  Thus  each 
particle  experiences  isothermal  deformations.  The 
relations  of  section  5.1  apply  to  each  point  but  the 
function  i|i  will  be  different  for  each  point  because 
the  temperature  is  different. 

It  was  expected  that  the  stress  at  the  inner 
boundary  due  to  a  step  change  in  boundary  pressure 
would  be  large  at  first,  then  diminish,  but  finally 
increasing  towards  the  initial  value:  As  shown  in 
Section  2.5,  the  initial  stress  should  be  the  same  as 
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for  an  elastic  material  and  therefore,  neglecting 
the  effect  of  temperature  on  the  value  of  G(o)  and  K(o), 
the  stress  should  he  the  same  as  in  the  uniform  tem¬ 
perature  test.  The  hotter  points  on  the  inner  boundary 
then  begin  to  relax.  The  modulus  G(t)  becomes  much 
smaller  at  the  hotter  points.  Consequently,  they 
"carry  a  smaller  portion  of  the  load",  and  the  stress 
will  decrease  at  the  hotter  points.  As  described  in 
Section  2.6,  the  stress  at  long  times  is  the  same  as 
for  an  elastic  material;  therefore,  neglecting  the 
effect  of  temperature  on  G(<»)  and  K(<»),  the  stress 
should  approach  that  observed  in  the  uniform  temperature 
test . 

The  fringe  order  actually  observed  at  a  point  is 
shown  in  Fig.  26.  At  a  point  on  the  boundary  the 
principal  axis  of  stress  is  normal  to  the  boundary  at 
all  times.  The  only  non- zero  stress  is  tangent  to  the 
boundary.  Choose  the  coordinate  axes  along  the  prin¬ 
cipal  axes  of  stress.  Then  (23)  of  section  5.1  applies: 

o  =  $ ( t —  t )  n( t )  dT  (2) 

o 

The  stress  calculated  using  the  observed  fringe  order 
and  the  function  §(t)  calculated  from  the  optical 
creep  function  i|i(t)  for  the  particle  temperature,  is 
shown  in  Fig.  27.  It  has  the  expected  character.  The 
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exact  analytical  solution  is  unknown  and  so  no  com¬ 
parison  is  possible. 
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(a)  loading  frame 

(b)  light  source,  filtered  Hg  arc 

(c)  Polarizer  and  l/4  \  plate 

(d)  tensile  specimen 

(e)  analyzer  and  1/4  X  plate 

(f)  photomultiplier  tube 

(g)  Brush  recorder 

(h)  thermoelectric  coolers 

(i)  power  supply  and  control  for 
coolers 

Figure  20  Tensile  calibration  bench. 
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Figure  3.  Optical  creep  compliance  Mix 
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Figure  6.  Poisson  ratio  in  creep,  Mix 
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igure  7.  Shear  creep  compliance  (calculated;,  Mix 


B(t)  =  3D(t)[l-2Z/(t)] 
- Z/(t)  =  Z/(oo) 
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al  expansion  determination. 


Figure  13.  Shear  model. 
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(a)  light  source,  filtered  Hg  arc 
(t>)  Polaroid,  rotating 

(c)  loading  frame 

(d)  load  cell 

(e)  model 

(f)  timer  and  angle  indicator 

(g)  Polaroid,  rotating 

(h)  camera 

(i)  rotating  Polaroid  control  and  power  supply 

(j)  drive  solenoids 

(k)  solenoid  valves  for  pressure  loading 


Figure  14.  Photovisocleastic  "bench. 
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SHEAR  STRESS  ALONG 
CENTERLINE  ( lb/in*) 


TEST  MODEL 


Figure  15.  Shear  stress  distribution. 
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Figure  16.  Fringe  order  in  shear  test 
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(a)  nichrome  heater 
(h)  latex  diaphragm 
(c)  guards,  (Polaroid  laminated  to 
Plexiglas) 


Figure  20.  Pressure  loading  jig  (model  removed) 


Figure  21.  Pressure  calibration  specimen  with  circula 
hole  (star  cutouts  in  Polaroid). 


Figure  22. 


(a)  Potentiometer 

(b)  Thermocouple  probe 

(c)  measuring  stations 

(d)  heater 

(e)  pressure  loading  jig 


Apparatus  for  determining  temperature  distribution. 
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(a) 

light  source,  filtered 

Hg 

(b) 

constant  pressure  air  s 

upply 

(c) 

loading  jig  (Polaroids 

incorporated ) 

(d) 

camera 

( e) 

heater  Variac 

Figure  23.  Horizontal  polariscope  for  heated  model  studies. 
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Figure  24.  Isothermal  pressurized  grain  model 
6300  seconds  after  loading. 
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Fringe  order  at  uniform  temperature 
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Figure  27.  Stress  at  star  root  (point  A),  heated 
model . 
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